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ABSTRACT 

We present the first ~7.5'xll.5' velocity-resolved (~0.2 kms -1 ) map of the [Cn]158/im line toward 
the Orion molecular cloud 1 (OMC 1) taken with the Herschel /HIEI instrument. In combination with 
far-infrared (FIR) photometric images and velocity-resolved maps of the H41 a hydrogen recombination 
and CO J- 2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinemat¬ 
ics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The 
main contribution to the [C n] luminosity (-85 %) is from the extended, FUV-illuminated face of the cloud 
(Go>500, ^h>5x10 3 cm -3 ) and from dense PDRs (Go>10 4 , tth^IO 5 cm -3 ) at the interface between OMC 1 and 
the Hii region surrounding the Trapezium cluster. Around -15 % of the [C ii] emission arises from a different 
gas component without CO counterpart. The [Cn] excitation, PDR gas turbulence, line opacity (from [ 13 Cn]) 
and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps 
of the L[C ii]/L F ir and L F iR/M Ga s ratios and show that L[C ii]/Lfir decreases from the extended cloud compo¬ 
nent (~10“ 2 -10" 3 ) to the more opaque star-forming cores (~10 _3 -10 -4 ). The lowest values are reminiscent of 
the “[C n] deficit” seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation 
analysis shows that the decreasing L[C ii]/Lfir ratio correlates better with the column density of dust through 
the molecular cloud than with Lpm/^Gas- We conclude that the [C n] emitting column relative to the total dust 
column along each line of sight is responsible for the observed L[C d]/L F ir variations through the cloud. 

Subject headings: Hu regions — galaxies: ISM — infrared: galaxies — ISM: clouds 


1. INTRODUCTION 

The 2 P 3 / 2 ~ 2 P \/2 fine structure line (AE/k^9l K) of ionized 
carbon, [Cn] 158yum, is one of the most important cooling 
lines of the cold neutral medium ( [Dalgarno & McCray|1972 ) 
and among the brightest lines in photocfTssociation regions 
(PDRs, e.g., lHollenbach & Tielens|1999| . Early observations 
showed that the line is very luminous ([Russell et al.|1980| ), 
carrying from - 0.1 to 1 % of the total fa r-IR (FIR) luminos¬ 
ity of galaxies (Craw ford et al. 1985). The neutral car¬ 
bon atom has an ionization potential of 11.3 eV, so that the 
ion C + traces the H + /H/H 2 transition, the critical conver¬ 
sion from atomic to molecular ISM. As a consequence, [C n] 
emission/absorption from different ISM phases is expected. 


Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with important 
participation from NASA. 

** Uses observations obtained with the IRAM-30m telescope. IRAM is 
supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain). 


Recent velocity-resolved line surveys show that the average 
contribution to the [Cn] emission in the Milky Way is the 
H 2 gas illuminated by far-UV (FUV; £^6-13.6 eV) photons 
(55-75 %), the cold atomic gas (20-25 %) and the Hu ionized 
gas (5-20 %) ( [Pineda et al.|20l3|[20T4| ). 


Star-forming complexes hosting massive OB stars are irra¬ 
diated by strong FUV fields and most of their [C n] emission 
is likely to arise from dense PDRs at the interface between 


Hu regions and their parental molecular cloud (e.g., Stacey 
et al. 1993| ). Therefore, the [C n] 158 jum line traces the FUV 


mation rate (SFR) (e.g 

., Stacey et al. 2010; Pineda et al. 2014;; 

Herrera-Camus et al. ! 

2015 ). As in the Milky Way, the [Cn] 


lated to the star formation activity ([Rodriguez-Fernandez et al. 
2006[ [Malhotra et al.|200T] Kapala et al.[20l5| ). In addition, 


an undetermined fraction of the [C n] arises from the diffuse 
neutral gas ( [Madden et al. 1993]) and from the H n ionized gas, 
the latter varying from ^5 to 50%, depending on the ion izing 
radiation strength and electron density (e.g., Abel|2006a| ). 
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In normal and starburst galaxies, the [Cn]-to-FIR lumi¬ 
nosity ratio (L[C ii]/L fir ) is in the 10 -2 to 10 -3 range. 
Low metallicity dwarf galaxi es sh ow the highest ratios 


(>10 , e.g., Mad den et al.||1997[ ), whereas observations 
with ISO showed that local ultra-luminous galaxies (ULIRGs, 
L F ir> 10 12 L o ) display a deficiency, L[Cn]/L FIR ^_lO~ 4 , that 
is not straightforward to interpret (Malhotr a et al.| [1997 ; 
Luh man et al.| [1998). This fractional luminosity de- 


ficiency has been confirmed with new observations of 
ULIRGs, hosting intense star formati on and often associ¬ 
ated with galaxy mergers (e.g., Gracia-Carpio et aT7||2011| 
|Gonzalez-Alfonso et al.[ |2015| ). This deficit, however, 
does not necessarily hold in the early Universe at high-z 
(e.g., |Stacey et al.|2010[|Brisbm et al.|2015| ). 

Interestingly, while it is difficult from the ground to detect 
the [Cn] emission in the local Universe (balloon, airborne, 
or space telescopes are needed to overcome the low atmo¬ 
spheric transmission at 158 /i m), the line does become acces¬ 
sible to ground-based submm observatories such as ALMA at 
redshifts z > L Indee d, the [C n] line has been detected toward 
young z> 5 galaxies (Riech ers et al.|2014 ; Capak e t al.|2015"] ). 

Most of the [Cn] emission arising from PDR gas in the 
Milky Way is associated with mode st average F UV radiation 
fields, between Go-1-20 ( Pineda et al.||2013] ) and G 0 -100 
( Cubick et al.||2008] ), with Go the mean interstellar FUV 
field in Habing units (Go=l is equal to 1.6xl0~ 3 erg cm -2 s -1 , 
|Habing| 19681 ). These radiation fields are likely not rep¬ 
resentative of the efficient massive star forming modes ex¬ 
pected in ULIRGs showing L[Cii]/L F ir deficits. Galactic 
high-mass star-forming regions, however, offer a more in¬ 
teresting template in which to investigate the [C n] emission 
in the more extreme irradiation conditions prevailing these 
active galaxies. The Orion molecular cloud 1 (OMC 1), in 
the Orion A c omp lex, lies behind the Orion Nebula cluster 
(e.g., [O’Dell|2001| ) and is the nearest (-414 pc), and probably 
most studied high-mass star-forming region of the Milky Way 
(e.g., |Genzel & Stutzki|l989| . 

OMC 1 is directly exposed to the UV radiation emitted by 
young massive stars in the Trapezium cluster, located ~0.3 pc 
in front of the molecular cloud (e.g., Bally|200>8 ). A region of 
several parsecs in size is exposed to FUV fluxe s much higher 
than the average Gp in the Milky Way (e.g., |Stacey et al. 


|1993[|Luhman et al.|1994[[Rodriguez-Franco et al.|1998| ). The 

Trapezium stars have created a blister Hu region that on the 

far side is bounded by the molecular cloud. The most famous 
Hn/OMC 1 interface is the Orion Bar PDR (the “Bright Bar” 
seen in the visible-light images), a protrusion of OMC 1 in 
which the neutral cloud acquires a nearly edge-on geometry 
and thus the optically thin PDR emission is li mb-brightened 
( [Hogerheijde et al. 1995 , Cuadrado et al.|2015| and references 
therein). On the near side, toward the observer, several com¬ 
ponents of neutral atomic gas rev ealed by H i absorption co ver 
the nebula, the Orion’s Veil (e.g., van der Werf et al.|2013 ). 

Intermediate- and high-mass star formation is cur¬ 
rently taking place in two locations: the Becklin- 
Neugebauer/Kleinmann-Low (BN/KL) region within OMC 1 
(Blak e et al.|1987[|Tercero et al.|2010| ), and Orion South (S), 

3 In this work we refer to the luminosity ( L ) computed in power units 
(ergs -1 or L©). The conversion from velocity-integrated line intensities 
f A7b dv = W(Kkms -1 ), with ATb the continuum-subtracted brightness 
temperature (see Section Lf, to line surface brightness (/), is I=2k Wv 3 /c 3 . 
For the [Cn] 158/im ana CO/=2-l lines respectively, the conversion is 
/(ergs -1 cm -2 sr -1 )=7.0xl0 -6 W (Kkms -1 ) and 1.3xlO -8 W (Kkms -1 ). 


an isolated molecular core located within the ionized gas cav¬ 
ity in front of OMC 1 (e.g., O’Dell et al.| 2009| ). There¬ 
fore, Orion is a unique region to study the spatial distribu¬ 
tion of [Cn] and related quantities (L[Cn]/L FIR , L[Cii]/Lco> 
L F iR/M Ga s ? etc.) in well characterized environments (dense 
PDRs, shocks, Hu regions, and lower density quiescent gas). 

The [C n] 158 /im line has been previously mapped toward 
OMC 1 at low angular (-55") and spectral (-67 km s -1 ) reso¬ 


lution with the Kuiper Airborne Observatory (KAO, St acey et 
|al.|1993[ [Herrman n et al.|1997]) and also w ith balloon-borne 
FIR telescopes (e.g.,|Mookerjea et al. 2003). Pointed obser¬ 
vations towardtheTr^SmmcIusterregion were also carried 
out by Boreiko & Betz|( |199~6| with a heterodyne receiver on 
board KAO. These observations employed a velocity resolu¬ 
tion of 0.5 kms' 1 and allowed the detection of two of the 
three [ 13 Cn] hyperfine line c omponents. The H IFI instru- 
ment ( de Graauw et al. 12010] ) on board Herschel ( [Pilbratt et| 


al.|20io ) allowed us to spectrally resolve (-0.2 km s 


reso¬ 
lution) ~and to map the [Cn] 158 /am line with unprecedented 
sensitivity and angular resolution (-11.4"). Velocity-resolved 
spectroscopy offers the best way to characterise the origin of 
all possible [Cn] emission components and their kinematics. 
In addition, it allows one to map [ 13 C n] and estimate the [C n] 
line o pacity, a critical parameter to determine C + column den 


sities ([Stacey et al. |1 99 lb[ [B oreiko & Betz|1996| [Goldsmith 


et al.|2012HQssenkopf et aLpOLpT 

In this work, we present initial results from a Herschel 
Open-Time program (OTl_jgoicoec_4) devoted to the spatial 
and velocity structure of the ionized and the warm molecular 
gas in OMC 1. In particular we present a large-scale velocity- 
resolved map of the [Cn] 158yum line, complemented by a 
map of the H41 a recombination line, and existing CO 7=2-1 
and FIR photometric observations. The paper is organized 
as follows. In Section [2] we present the data set. The [Cn] 
spatial distribution and gas kinematics is studied in Section]!] 
In Section [4] we analyze the [C n] 158 jum non-LTE excitation 
and derive C + column densities. A spatial correlation analy¬ 
sis of different quantities often used in the extragalactic dis¬ 
cussion is also carried out. Finally, our results are discussed 
in the context of the [C n], CO and FIR diagnostic power and 
their link with the extragalactic emission. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. [Cn], H41a and CO 2-1 Heterodyne Observations 

The [Cn] line at 1900.537 GHz ( Cooksy et al.||1986| ) was 
mapped over a total area of ~7.5 , xll.5' (~0.9pc xL4pc) in 
2012 August, using the HIFI instrument on board Herschel. 
A pair of Hot Electron Bolometer mixers was employed 
to observe the [Cn] line in HIFI band 7b in two orthogo¬ 
nal polarizations. We used the Wide Band Acousto-Optical 
Spectrometer providing a spectral resolution of 1.1 MHz and 
a bandwidth of 2.4 GHz (-0.2 km s _l and -380 km s _l at 
1900.5 GHz, respectively). 

Three ~7.5'x4 / On-The-Fly (OTF) submaps of 3h each 
were needed to cover the region with a 10" sampling 
(ObsIDs: 1342250415, 1342250414 and 1342250412). The 
half-power beam width (HPBW) of HIFI at these frequen¬ 
cies is -11.4", i.e., the map is not fully_sampled. Because 
of the large scale [Cn] emission ( [Stacey et al.||1993| and 
risk of contaminating emission in theTeference position, the 
OTF ma ps w ere performed using the Load-Chop referencing 
scheme (Roelfsema et al. 2012). In this mode, an internal 
chopper alternates between the sky and an internal cold load 
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Figure 1. (Left): Composite image with the H41a (green), [C n] 158 //m (blue) and CO 2-1 (red) integrated line intensities. The positions of the main sources in 
OMC 1 are shown. (Right): Herschel/YHFI map of the total [Cii] 158 gm integrated line intensity 3 (in Kkms -1 ) from vlsr=-30 to +30 kms -1 . 


Table 1 

Spectroscopic parameters of the lines discussed in this work 


Species 

Transition 

Frequency 

(GHz) 

EJk 

(K) 

Aii 

(s' 1 ) 

12 c + 

2 P3/2~ 2 Pl/2 

1900.5369“ 

91.21 

2.32(-6) 

13 C + 

2 Ps/2~ 2 Pl/2 
F= 1-0 

lgOO^OO"’'’ 

91.23 

1.55(-6) 


F= 2-1 

1900.4661“' ( ’ 

91.25 

2.32(-6) 


F= 1-1 

l900A360 a ’ b 

91.23 

7.73(-7) 

CO 

/=2-1 

230.5379 

16.60 

6.92(-7) 

H41 a 

n=42-41 

92.0344“ 

89.50 

4.85(+l) 


Note. —r a From Cooksy et al. 

(1986 

i. ^SeelOssenkopf et al. ( 2013) 

corrected relative line-strengths. L J 

From 

.dlley & Palmer (19681. 


as the telescope slews over a scan leg. Although a nearby 
reference position (at a radial offset of - 9' from the map 
centre) was also observed, we reprocessed the data without 
subtracting this reference emission from the on-source data, 
therefore being insensitive to line contamination from self¬ 
chopping. The resulting data are mostly affected by so-called 
Electrical Standing Waves, that were removed in HIP^] us¬ 
ing the doHebCorrection task ( Kester et al. 2014| ). The resid¬ 
ual optical standing waves resulting from resonances against 
the internal hot load (-100 MHz period) and the diplexer op¬ 
tics (-625 MHz) were treated using the FitHifiFringe task in 
HIPE. The three sub-maps were then compared pair by pair 
in a region of -7.5'xO.6' where they overlap on the sky. An¬ 
tenna temperatures (7^) were extracted for each overlapping 
pixel and [C n] velocity channel. Scatter plots were formed 
to derive the inter-calibration factor between the respective 
coverage in each channel. We derived linear slopes deviat¬ 
ing from unity by -5% between pairs of sub-maps. Since we 
cannot know a priori which map corresponds to the true in¬ 
tensity, we re-adjusted the fluxes by ±2.5% depending on the 
sub-map, i.e. re-calibrated with respect to the mean intensities 
in the overlapping regions. 


Through this work we use the main beam temperature scale 
(T mh in K) as opposed to T* h . Given the latest HIFI effi¬ 
ciencies, both scale^are related by r mb -1.64 T* A . For semi- 
extended (not infinite but larger than the main beam width of 
Herschel/ HIFI) emission sources of uniform brightness tem¬ 
perature (7b), the main beam temperature is the most appro¬ 
priate intensity scale (T mb ^r b ). The achieved rms noise was 
typically -1.5 K (l<x) per 0.2kms -1 resolution channel and 
position. T his rms is only a factor of -2.5 worse than that 

in their deep pointed ob¬ 
it. Figure [I] (right) shows 


achieved by Boreiko & Betz 
servations after an integration^ 

the baseline-subtracted [C n] integrated line intensity map. 

The H41 a millimeter recombination line at 92.034 GHz 
was mapped with the IRAM-30m telescope at Pico Veleta 
(Spain) using the E090 receiver and the 200 kHz FFTS back¬ 
end in 2013 November. The HPBW at this frequency is -27" 
and the typical rms noise in the map is low, -0.05 K per 
0.65 kms -1 resolution channel. In this work we also make 
use of part of the high-sensitivity CO 2-1 large-scale map 
(1x0.8°) at 230.538 GHz obtained by|Berne et al.|(|2014| with 
the multi-beam receiver HERA, also at the IRAM-30m tele¬ 
scope. The spectral and angular resolutions are 320 kHz and 
-11", respectively. The achieved rms noise is -0.2 K per 
0.4kms -1 channel (see Berne et al. 2014). Both H41 a and 
CO 2-1 OTF maps are fully sampled. 

Herschel /HIFI and IRAM-30m data were processed with 
the GIFDAS/CFASS software. Figures [lj[6j and [7] show the 
[Cii] 158//m, H41 a and CO 2-1 emission at their native an¬ 
gular resolutions but resampled to a common velocity grid 
of 0.65kms -1 . Offsets of selected positions are given with 
respect to the center of the [Cii] map: o^ooo- 5"35 m 17.0 5 , 
£ 2000 • -5°22 , 33.7". Table [I] summarizes the main spectro¬ 
scopic parameters of the lines discussed in this work. To 
match the FIR photometric observations and carry out a com¬ 
bined analysis, the line maps were also gridded and convolved 
to an uniform angular resolution of 25" (-0.05 pc). 


4 HIPE is a joint development by the Herschel Science Ground Segment 

Consortium, consisting of ESA, the NASA Herschel Science Center, and the 5 See: http: //herschel. esac. esa. int/twiki/pub/Public/ 

HIFI, PACS, and SPIRE consortia. HifiCalibrationWeb/HifiBeamReleaseNote_Sep2014 .pdf 
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[CII]158um (total) [CII]158um [-2,5] km/s 



Figure 2. Color-composite visible image of the Orion Nebula (M42) taken with the ACS instrument on board the Hubble Space Telescope (HST). The assigned 
colors are: red and orange for [Sn] and Her, blue for B, and green for V filters respectively (Robberto et al. pot's) . (Left): Green contours show the integrated 
intensity (W) of the H41er recombination line (from 3 to 21 K km s -1 in steps of 2 K km s -1 ). The white dotted 200" (0.4 pc) radius circle is centered on 9 l Ori C, 
the brightest star of the Trapezium. We define the “extended cloud face” component to lie outside this circle. Middle : Black contours show the total [C ii] 158 pm 
integrated intensity (vlsr=-30 to +30 kms -1 with W from 0 to 1300 Kkms -1 in steps of 150 Kkms -1 ). Right : [Cii] “blue-shifted streamer” (emission from 
vlsr- - 2 to +5 km -1 ) with W from 0 to 60 Kkms 1 (cyan) and from 80 to 140 Kkms 1 (blue), all in steps of 20 Kkms -1 . This emission coincides with the 
Northern dark lane and with the Dark bay seen in visible-light images of the Nebula. The Orion hot core in the BN/KL region and Orion S positions are indicated 
by red and blue stars respectively. The positions where line spectra have been extracted (shown in Figure[3j are indicated by cyan squares. 


2.2. Photometric Observations and SED Fits 

In order to determine the FIR luminosity and the dust con¬ 
tinuum opacity toward each position of the [C n] map, we also 
make use of fully sampled calibrated images of OMC 1 taken 
with the Herschel /SPIRE an d PACSQ cameras ( Griffin et al. 
2010} [Poglitsch et aL]|2010| ) at 70, 160, 250, 350/rm photo¬ 
metricTan3s"(Ifonrprograms GT2_pandre and OTl_nbillot), 
and with JCMT/SCUBA2 at 850yum (project JCMTCAL; 
[Holland et al.|2013) . 

We convolved all the photometric images to a 25" resolu¬ 
tion (that of the 350/mi image) and fitted the resulting spec¬ 
tral energy distribution (SED) per 12.5"-pixel with a modified 
black body at an effective dust temperature (T d ); 

Fa = BffT d ) • (1 - e~ T ^) • Q (1) 


where Q is the solid angle subtended by each pixel. Fit exam¬ 
ples are shown in the Appendix. The dust continuum opacity 
is parametrized as r d j = T d ^(\6QI Xf, with r d ,i6o the dust 
opacity close to the [Cn] 158 pm line, and /3 a grain emis- 
sivity index. The resulting L FIR and r dj i6o maps, where we 
have allowed to be a free parameter of the fits, are shown 
in Figures[8ja) and (f). The mean value is 2.4+0.5 (l<x) 
but we note that the derived FIR luminosities and dust opac¬ 
ities, the relevant quantities for our analysis, do not change 
significantly by fixing /3 to a constant value of 2. We also es¬ 
timated the pixel-averaged gas column density and mass. In 
particular, we used the SPIRE 250 pm image assuming that it 
corresponds to optically thin dust emission and computed the 
equivalent gas mass per pixel as: 


^Gas = 


Rg&S 250 F> 2 
*250^250 (Td) 


( 2 ) 


( [Hildebrand] 1 1983 1 ) , where S 250 is the 250/rm flux density 
within Q, A gd is the gas-to-dust mass ratio and ^250 is the dust 


6 Lfir from 40 to 500 pm ( Sanders & Mirabel 1996} and using D= 414 pc. 

7 The central region of the photometric image~(BN/KL outflows and IRc 
sources) contain saturated data points that are not used in our analysis. 


absorption cross-section divided by the dust mass. We adopt 
^gd/k25o=25 grem -2 , computed from grain properties taken 
from Draine’s tabulations ( |Li & Draine|2001[ and references 
therein) for a grain size distribution that leads to a A v /Nu ratio 
compatible with observations of Orion (see Section |L2|). 


3. RESULTS 

3.1. Global Morphology of the [C u] 158pm Emission 

Figure [T] (left) shows a composite image of OMC 1 in the 
integratealine emission of CO 2-1 (red), [Cn] 158 pm (blue) 
and H41a (green). This image shows the H 11 /PDR/molecular 
gas stratification expected in molecular clouds irradiated by 
strong UV fields from nearby massive stars. 

The H41 a recombination line traces dense ionized gas from 
the blister Hu region surrounding the Trapezium. The H41 a 
line peaks near the OB stellar cluster and is also bright ad¬ 
jacent to the Orion Bar PDR. Despite the different angular 
resolutions, the H41 a emission morphology is similar to the 
brightest regions seen in the visible-light H a images of the 
Orion Nebula (M42) (Figure |2| /g/f; [Robberto et al.|2013| ). 

Figure [3] shows [Cn], CO 2-1 and H41<+ spectra extracted 
from representative lines of sight toward OMC 1, namely the 
northern streamer, BN/KL, the Trapezium, the Orion Bar, the 
East PDR, and the extended cloud face in the south. OMC 1 
has an overall face-on geometry, with the molecular cloud lo¬ 
cated ~0.3-0.6 pc behind the ionizin g stars (e.g., Hogerheijde 
|et al. 1 1995} | van der Werf et al.|2013| ). Therefore, the observed 
— 10-15 km s -1 blue-shift of the H41 a line relative to the [C 11 ] 
and CO lines (Figure[3]) is consistent with streamers of ionized 
gas escaping fromthenigh-pressu reH 11 /OMC 1 interface and 
moving toward the observer (e.g., |Genzel & Stutzki|1989| ). 

The regions of brightest [C 11 ] 158 pm emission delineate 
the far edges of the H 11 blister and follow a roughly spherical 
shell distribution (as projected in the plane of the sky) with 
emission peaks near the Trapezium, at the east of BN/KL (that 
we refer to as the “East PDR”) and along the Orion Bar PDR 
(Figure [T] right). These “interface” regions show very bright 
[C 11 ] emission, with main beam temperature peaks of >150 K 
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Figure 3. [Cn] 158/mi, CO 2-1 and H41 a line profiles toward representa¬ 
tive positions in Orion. Offsets in arcsec are given with respect to the HIFI 
map center at a 2 00(F 5*35 OT 17.0 s , 62000 ' -5°22'33.7". The last panel shows 
the average spectrum over the region. The inset panels show a zoom into 
the spectra. The magenta arrows mark the expected positions of the [ 13 C 11 ] 
F=l-0, 2-1 and 1-1 hyperfine components (from left to right). The [Cn] and 
CO spectra were extracted from maps convolved to 25" resolution. 



(or a surface brightness of >5xl0 -3 ergs -1 cm -2 sr -1 for the 
average Av^5kms -1 line-width, see Figure pj). In addition, 
a more extended component of fainter [C n] lT8 pm emission 
exists. In particular, the average line peak temperature in the 
map is ~70 K. This reflects the presence of a lower density 
and cooler component, the extended FUV-illuminated “face” 
of OMC 1 previous ly detected in [C n] at lower angular and 
spectral resolutions jStacey et al .[|1993 ; Herrmann et al. 1991) . 
This extended component is also seen in the wide-field images 
of M42 taken in the emission from other low ionization po¬ 
tential (<13.6 eV) atoms (see the reddish nebular emission in 
Figure[2]taken with the [S n] filter in the visible). It is also seen 
in the near-IR vibrationally excited H 2 emission^attributed to 
large-scale FUV fluorescence ( [Luhman et al.|1994 ). 
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Figure 4. Spitzer/IRAC 8 pm filter image (dominated by PAH emission) 
toward OMC 1. Contours show the total [Cn] 158/mi integrated intensity 
(from 0 to 1300 Kkms -1 in s teps of 150 Kkms -1 ). The different circles and 
regions are defined in Section |T2| The positions where line spectra have been 
extracted (shown in Figure[3j are indicated by squares. 


Most of the [C ii] (and obviously CO) emission is expected 
to arise from the neutral cloud. This can be readily seen 
from the similar [Cn] 158 yum and CO 2-1 line profiles and 
also similar line peak velocities (Figure [3]). At the small spa¬ 
tial scales revealed by our observations (-5000 AU), how¬ 
ever, both the integrated intensity map (Figure [I] left) and 
the velocity channel maps (Figure [7]) show that the [Cn] and 
CO 2-1 spatial distribution does not exactly match. This is 
more obvious toward (nearly) edge-on PDRs, like the Orion 
Bar, because the irradiated cloud edge (emitting [C n]) can 
be resolved from the more shielded molecular gas, where 
CO reaches its abundance peak (A v >3 mag for dense gas and 
strong FUV fields, see the representative model in Fig. [5]). 

The integrated intensity maps also reveal regions of moder¬ 
ately bright [Cn] 158yum emission, but very weak CO emis¬ 
sion counterpart (northeastern areas). In addition, they show 
regions, such as the shocked gas toward H 2 Peaks 1 and 2 in 
Orion BN/KL outflows, wh ere the CO line lumino sities are 
exceptionally high (see also Goicoechea et al .j2015| . 

Finally, Figure [4] shows the good spatial correspondence 
between the [Cn] 158 yum integrated i ntensity (contours) and 
the Spitzer/IRAC 8/rm emission (e.g., |Megeath et aL|2012| ). 
This image is dominated by fluorescent IR emission from 


PAHs (af ter absorption of UV photons, e.g., Allamand ola et 
al. 1989) and by very small grain emission, with minor con 


tribution from [Ar hi] toward the H n region and from H 2 lines 
toward particula r pos itions of enhanced shocked gas emission 
( [Rosenthal et al.|2000[ ). 


3.2. Template Regions and Representative PDR Model 

Based on our knowledge of Orion (e.g.,|Genzel & Stutzki 
1989[ |Stacey et~aLl 1993 Bally||20Q8] ), and in order to ease 
the interpretation of the \C nl and relate 


i interpretation of the [C n] and related maps, we introduce 
three different template regions (or sightlines). We first con¬ 
sider a 200" (0.4 pc) radius circle centered on 0 l Ori C, the 
brightest star of the Trapezium (e.g., the white dotted cir¬ 
cle shown in Figure |4]). We define the “extended cloud face” 
component for positions outside this circle. Secondly, a small 
R<0.1pc circle around the Trapezium defines lines of sight 
that go through the “H n region” component around the H41 a 
emission peak. Finally, we consider a spherical shell, or annu- 
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Figure 5. H/H 2 and C + /C/CO transitions predicted by the Meudon PDR code 
for physical conditions and grain properties representative of OMC 1 near 
the Trapezium cluster (Go=2xl0 4 andnH=n(H)+2n(H2)=2xl0 5 cm -3 ). The 
equivalent length of the Ay= 10 mag slab is ^10,000 AU (^24" at the distance 
to OMC 1). The dotted C + and CO abundance profiles in the uppe r pan el are 
obtained using standard, diffuse ISM, grain properties (see Section|+2). 


lus, delineated by the R=0.16pc and R= 0.24 pc circles. This 
region represents “PDR interfaces” between OMC 1 and the 
Hu region. These PDRs (specially the Bar) show a more 
edge-on configuration and are illuminated by the strong, al¬ 
most unattenuated, UV field from the Trapezium stars. 

To guide our interpretation, Figure [5] shows the predicted 
structure of the H/H 2 and C + /C/CO transitions as a function 
of visual extinction (Ay) for physical conditions (Go=2xl0 4 
and /?h=2x10 5 cm -3 ) representative of the Hn/OMCl in¬ 
terface near the Trapezium (see Section [4]). This homo¬ 
geneous, stationary slab cloud model was computed with 
the Meudon PDR code (Le Pe tit et ah||2006| . We adopted 
[C/H] = 1.4x 10 -4 (from HST absorption obser vations tow ards 
the star HD 37021 in the Trapezium, [Sofia et al.||2004| ) and 
a flat extinction curve compatible with observationsof Orion 
dLee|| 1968| |Cardelli et alJ[T989l |Allers et al.||2005l |Abel et| 
al. 2006b| ). In particular, an extinction to color index ra¬ 


tio (Ry = Av/Eb-v) of 5.5, and Ay/AH=3.5xlO _22 magcm 
were adopted. These values are different than the standard 
grain properties used in diffuse cloud models (Ry = 3.1 and 
A v /Nu=53xl0~ 22 mag cm -2 ) and are consistent with larger 
than standard size grains in Orion. Larger grains lead to an 
increased penetration of FUV photons, larger size of the C + 
zone (AAy(C + )>l mag), higher C + column densities, and shift 
the peak CO formation layer to larger A v deeper inside the 
cloud (see also |Goicoechea & Le Bourlot|2007| ). 

More in general, an ionization front can move into the 
molecular cloud. Ionized gas will rapidly photoevaporate 
from the cloud surface and the PDR structu re will not be 
in equilibrium (e.g., Bertoldi & Drain e|1996] ). Thus, dy¬ 
namical effects might alter the picture shown in Figure 
Non-stationary PDR models predict that a well define* 


1 


C + /CO transition layer always exists, although the layer is 
slightly shifted closer to the cloud surface (by only AAy^0.5 
for Go/ftH conditions appropriate to OMC 1) due to advection 
of CO from inside the cloud ( [Storzer & Hollenbach||1998] ). 
Since the predicted C + column density in these models does 
not vary much compared to that in stationary models, the gen¬ 
eral predictions shown in Figure[5]for C + should remain valid. 

3.3. Kinematics: [Cn] 158 pm Velocity Channel Maps 

Figures [6] and [7] show continuum-subtracted [Cn] 158 pm 
velocity cnannel maps compared to those of H41 a and 
CO 2-1 lines respectively. Both figures show the emission 
in 1 kms' 1 velocity channels from vlsr = -9 to +20 kms -1 
([Cn] 158 pm in colored maps and H41a or CO 2-1 lines in 
white contours). They reveal the intricate small spatial-scale 
kinematics of the ionized/PDR/molecular gas interfaces. 

The velocity-centroid of the molecular cloud emissioij^] 
traced by CO, liesatv L sR -+(8-10) kms -1 . The [Cn] 158 yum 
spectrum average over all mapped positions is characterized 
by a very bright component peaking at Vlsr - +9.5kms -1 , 
with a line-width of ~5 kms -1 (lowest panel in Figure [ 3 }, in 
agreement with CO velocity and line widths. 

This emission mostly arises from the H 11 /OMC 1 interfaces 
close to the Trapezium (inside the 200" radius circle) and also 
from the large-scale FUV-illuminated face of the molecular 
cloud (outside). Therefore, the main [Cn] spectral compo¬ 
nent that coincides with CO (vlsr ~5-17kms~ l ) is dominated 
by C + in dense PDR gas and represents ~85 % of the total 
[C 11 ] luminosity in the map. The emission is so bright that the 
[ 13 Cn] F= 2-1 line, the strongest [ 13 Cn] hyperfine component 
( Ossenkopf et al.||2013| ), is also detected toward the highest 
column density peaks (see the inset in Figure [3] and the chan¬ 
nel maps between +17 and +20 kms -1 in Figures [6] and Lj). 

Other specific regions such as the Orion Bar PDR show 
[Cn] line peak velocities at vlsr -+(10-12) kms -1 and thus 
have a slightly differentiated kinematics. In fact, [C 11 ] shows 
widespread emission in a broad range of velocities outside 
the main spectral component range. On the other hand, 
the CO 2-1 emission range is much more restricted, with 
the high velocity CO emission seen at VLSR<0kms -1 and 
vlsr>15 km s -1 being very localized, and mostly arising from 
the molecular outflows surrounding Orion BN/KL and S. 

Despite the relatively flat [Cn] intensity distribution sug¬ 
gested by previous lower resolution observations, the [Cn] 
channel maps show complicated structures and velocity pat¬ 
terns. This is true even away from the Trapezium, where fila¬ 
mentary structures and emission striations are seen at several 
velocity intervals. At the vlsr ~+8-10 kms -1 velocity of the 
neutral cloud, the H41 a emission delineates the far edges of 
the expanding H 11 blister in close interaction with the cloud. 
This produces a bright shell-shaped [C 11 ] emission from com¬ 
pressed PDRs immediately beyond the ionization fronts. This 
is the signature of the UV radiation eating into the cloud. 

Also at vlsr ~+8-10 km s -1 , but away from the H 11 region, 
[C 11 ] shows narrow striations or sheet-like structures that are 
roughly perpendicular to the main molecular ridge (north-to- 
south) and do not show much CO emission. As suggested by 
the overlap with the visible-light images of M42 (Figure [2]), 
these photoevaporative [C 11 ] structures demonstrate that the 
entire field and velocity range is permeated by FUV photons. 

The most negative velocities show H41 a emission around 


; The Orion A molecular complex to which OMC 1 belong s show s a ve¬ 


locity gradient of 0.7 km s 1 pc 1 over ~25 pc (e.g., Bally et al. 


19871. 
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Figure 6. Velocity channel maps in steps of 1 km s -1 from vlsr = -9 to +20 km s -1 ([C n] in colors and H41 a line in contours). The [C ii] main beam temperature 
scale in bins of 1 kms -1 is given in the color wedge (in some velocity intervals the [Cii] intensity per channel is multiplied by a factor that is specified in each 
panel). Offsets in arcsec are given with respect to the HIFI map center. The H41cr contours go from 0.2 to 1.8 Kkms -1 in steps of 0.3 Kkms -1 . The [Cn] 
emission seen at vlsr = +20kms -1 (last panel) is dominated by the [ 13 Cn] F= 2-1 line peak. Note that the first two upper-left panels show total integrated line 
intensity maps for H41a (from 0 to 22 Kkms -1 ) and [Cii] 158 jum (from 0 to 1300 Kkms -1 ). In the total [Cii] intensity map, H41a integrated line intensities 
are shown as white contours (from 1 to 21 Kkms -1 in steps of 4 Kkms -1 ). The Orion hot core, Orion S condensation and Trapezium stars are indicated by red, 
yellow and blue stars respectively. The positions where line spectra have been extracted (shown in Figure[3j are indicated by white squares. 


the Trapezium. This is probably the edge of the Hu blister 
that is expanding toward the observer. The analogous [C n] 
emission has some spatial correspondence but it also shows 
specific blue-shifted structures that do not follow the H41 a 
emission and that must lie in the line of sight toward OMC 1. 

3.3.1. C + Gas with no CO Counterpart 

In certain velocity ranges the observed [C n] emission does 
not have any CO emission counterpart. The most outstanding 
structure with no, or very little, CO emission is a blue-shifted 
[Cii] emission “streamer” (vlsr=-2 to +5 km -1 ) that goes 
north of Orion BN/KL and also a compact region that fills the 
cavity inside the bright [Cii] shell-like structure (Figure [2]). 
These structures are not coincident with H41 a emission trac¬ 
ing H ii gas, suggesting that they arise from warm neutral gas 
of low shielding (Ay^0-4mag), in which very little or no CO 
exists (due to the elevated FUV field in the region and thus 
rapid photodissociation). Interestingly, they do coincide with 
foreground depressions of the visible-light emission known as 
the Dark Bay and the Northern Dark Lane that separates M42 
and M43 nebulae (blue contours in Figure [2] right). 

The column density of material in these structures is suf¬ 
ficient to have significant hydrogen in molecular form (con¬ 


sistent with the detection of rad io OH lines in abso rption to 
ward the Dark Bay region, e.g., |Brogan et al.|2005|. Hence, 
they can be referred to CO-dark H 2 gas ( |Grenier et aL[ 2005 


Wol fire et al.||2010| ), whi ch can also be d etected b y dus t ex 


tinction of visible light ( [O’Dell & Yusef-Zadeh 2000| ). In 
[Cii] emission, they represent large-scale structures of C + 
gas surrounding OMC 1 and showing blue-shifted veloci¬ 
ties relative to the molecular cloud (dark blue contours in 
Figure [2] right). In terms of luminosity, the [Cii] emission 
from vlsr - -2 to +5 km s _1 , not related to OMC 1, accounts 


for ~12% of the total [Cii] luminosity. 

Surprisingly, the [C 11 ] line shows emission components at 
even more negative LSR velocities throughout the region, in 
particular in the Vlsr - -10 to -2kms -1 range, and a fainter 
blue emission wing at vlsr - -25 to -10kms -1 (see spectra 
in Figure [3]). Their combined contribution to the total [Cii] 
luminosity is however small, ~3%. 

Some of the structures showing [C 11 ] emission but no CO 
counterpart can be associated, both spatially and in veloc¬ 
ity space, with particular regions of H 1 absorption observed 
against the radio continuum. They are produced by the “Veil”, 
a collection of absorbing H 1 layers in the line of sight toward 
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Figure 7. Same as Figure[6]but the CO 2-1 line (in contours). The CO 2-1 contours go from 1.5 to 160 Kkms -1 in steps of 30 Kkms -1 . The high velocity CO 2-1 
emission (vlsr<0 s -1 and vlsr>15 s -1 ) mostly arises from molecular outflows around Orion BN/KL and Orion S protostellar sources. At other positions, the 
observed [C n] emission in these velocity ranges does not have CO emission counterpart (e.g., the blue-shifted [C ii] streamer at vlsr=-2 to +5 km -1 , or the [C n] 
emission in the southwest at vlsr ^ 12kms -1 beyond the Orion Bar). The first two upper-left panels show total integrated line intensity maps of CO 2-1 (from 
0 to 600 Kkms -1 ) and [Cii] 158//m (from 0 to 1300 Kkms -1 ). In the [Cii] total intensity map, CO 2-1 integrated line intensities are shown as white contours, 
from 200 to 500 Kkms -1 in steps of 100 Kkms -1 (extended emission) and from 800 to 1400 Kkms -1 in steps of 200 Kkms -1 (Orion BN/KL outflows). 


Orion with little H 2 , les s than one part in 10 6 ( Troland et al. 
1989; | Abel et ak1|2006b| ). For a detailed explanation of the 
different ~Hi ab sorbing velocity com ponents (A, B, C, D ...) 
see e.g., |van der Werf et ah] ( |2013| ). In particular, compo¬ 
nents A and B of the Hi absorption maps, produce strong 
Hi opacity in the Vlsr-0 to +6kms _1 range. Components A 
and B are interpreted as a foreground veil of nearly atomic 
gas in front of the Orion nebula. Thus the [Cii] emission 
in the vlsr-0 to +6kms _1 range, and outside the Dark Bay 
and North Dark lane regions, likely arises from these neutral 
atomic layers of the Veil (cyan contours in Figure [2] right). 

Component C of the H i absorption (from 
vlsr- - 4 to -2kms _1 ) is only found at the southwest¬ 
ern edge of OMC1 and we also detect [C n] emission 
at those velocities and positions. Finally, Hi absorp¬ 
tion components D1 and D2 near the Bar appear in the 
vlsr--20 to -lOkms -1 range and are interpreted as an 
expanding shell of atomic gas ( van_der Werf et al j20 13| ). 
The most blue-shifted [Cii] emission (vlsr<-10 knis -1 ) is 
detected toward the southern regions of the map. Hence, 
part of this [Cii] emission may arise from D1 and D2. In 
addition, a fraction of this blue-shifted [Cii] emission may 
arise from the foreground Hu region that lies between the 


Veil and the Trapezium, and that is not the main H n region 
interacting with OMC1. [Nn] T6585A emission line, and 
Pm T1335A and Sin /U190A absorption lines centered 
at vlsr<-13 kms' 1 arise from this foreground Hu region 
(e.g., Abel et al.|2006b| ). Unfortunately, the sensitivity of our 
[C n] observations is insufficient to map these faint features 
and to determine a definite association. 

H i emission tracing denser gas associated with the molec¬ 
ular cloud is detected at Vlsr> +10 km s _1 , outside the strong 
Hi absorption range of the Veil. Indeed, we find specific 
regions that show coincident [C n] and H i red-shifted emis¬ 
sion. The brightest one is the Orion Bar, where both [C n] and 
H i trace the most exposed PDR layers, in which H 2 is pho- 
todissociated. In addition, an elongated [C 11 ] and H 1 emis¬ 
sion feature without CO counterpart (vLSR^+12kms -1 ) ex¬ 
tends in the southwest direction beyond the Bar and along the 
Bar prolongation seen at visible wavelengths (Figure [2]). An¬ 
other remarkable emission feature is the high-velocity [C 11 ] 
and Hi line-wing detected toward Orion BN/KL outflows. Its 
presence confirms that the shocked material is illuminated by 
FUV radiation ( [Chen et al.|2014[ [Goicoechea et al.|2015] ). 

All in all, our observations demonstrate the link between 
[C 11 ] and both H 1 absorption (hydrogen mostly atomic gas) 
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and Hi emission (hydrogen predominantly molecular). In 
particular, - 15% of the total [Cn] luminosity does not have 
a CO counterpart. Most arises from cold Hi gas and from 
CO-dark H 2 gas (e.g., the Dark Bay), but also from Hu gas. 
The contribution from ionized gas is expected to decrease 
with increasing electron density (n e ) and effective temper¬ 
ature (T e ff) of the ionizing stars (as carbon becomes dou¬ 
bly ionized). Both quantities are high toward the Trapez - 
ium, L eff ^39,000 K for G l Ori C star ( Simon-Diaz et al. 2006), 
and < n Q >^(0.5-l )xl0 4 cm -3 in the surrounding Hu region 
(Zuckerman [1 97 3 ) . For these conditions, models predict that 
the contribution of [C 11 ] from ionized gas i s small, of th e order 
of <10% of the PDR plus Hu emission ( |Abel||20d6a| ). This 
agrees with our general conclusion that most of the observed 
[C 11 ] luminosity in Orion arises from dense PDR gas. 

3.3.2. [Cn] Line Broadening and C + Gas Turbulence 

The [Cn] 158yum line-width of the main spectral compo¬ 
nent varies between Avfhwm~3.5 km s -1 toward the extended 
“face” of OMC1, -4.0 kms' 1 toward the Orion Bar, and 
-5.5 kms -1 toward the Trapezium. The two latter value s are 
likely affected by opacity broadening (see Section 4 A for a 
determination of [Cn] opacities). In particular, for T[cii]-l-3, 
opaci ty broadening increases the line widths by -10-30% 
(e.g., [Phillips et al.||1979[ |Ossenkopf et al.||2013| ). We note 
that the line widths of the C91 a radio recombinati on line to - 
ward the Orion Bar are 2.0-2.5 km s -1 (Wyrowski et al. 1997). 

In the dense PDR gas traced by the [C nj main spectral com- 
ponent inside the 200" radius circle, the g as tem perature is 
T k ^300-1000K (see Figure [5] and Section |4 1.1 1 ), implying 
a thermal broadening of <x t h=(kr k /mc+) 1 -0.5-0.8 kms -1 . 
The nonthermal velocity dispersion <x tur b (from turbulence 
and macroscopic motions) is cr tmh =(cr 2 - <x 2 h ) 1/2 where 
<t=Avfhwm/2.355. For the above line-widths and gas 
temperatures, one obtains cr turb =1.2-1.7 kms -1 (after opac¬ 
ity broadening correction). This has to be compared 
with the isothermal sound speed in the warm PDR gas, 
cpDR=(£?kM) 1/2 -l.4-3.0 kms -1 (with m the mean mass per 
atom). In other words, the nonthermal velocity dispersion 
in the dense PDR gas is subsonic or only weakly supersonic 
(with Mach numbers <x tur b/cpDR of ^0.4-1.2). 

In the intermediate Vlsr=-2 to +5 km -1 components, the 
[Cn] line is narrower, for example Avfhwm~2. 5-3.0 kms -1 
toward the blue-shifted streamer that coincides with the 
North-Dark-Lane. On the other hand, the most blue- 
shifted [Cn] component with neither CO nor H41 a coun¬ 
terpart, for example the features at VLSR--7kms _1 toward 
(-150",-300") or at v LS r^-5 Ians' 1 toward (-37",-274"), 
show again broader line-widths of Avfhwm~4.7 kms -1 
(Figure]?]). Assuming optically thin [Cn] emission and 
L k <300 K, we find an increased <x tur b-2kms _1 (or 
cr turb/^Hi>2) in these CO-dark components surrounding 
OMC 1. These results suggest that the gas turbulence is more 
strongly dissipated in the dense H 11 /OMC 1 interfaces than in 
the surrounding, more turbulent, halo components. 

3.4. Surface Luminosity Maps 

Figure M shows maps of different luminosit^and luminos¬ 
ity ratios that are often discussed in the context of star forma¬ 
tion and galaxy properties (e.g., Carilli & Walter|2013[ ). 


3.4.1. Dust Opacity, Mass and Surface Density 

Figure [8ja) shows a map of Lfir per pixel obtained by 
integrating the dust SED 6 fits at each position. Despite 
the relatively small area covered (-1.25 pc 2 ), the FIR inten¬ 
sity exhibits variations of -3 orders of magnitude (excluding 
Orion BN/KL). Lfir is obviously high toward the high column 
denshy pl regions in the molecular ridge (Ah > 10 23 cm -2 ) 
and locally, it reaches even higher values toward the em¬ 
bedded star forming cores in Orion BN/KL. Figure [8T) 
shows a map of the dust opacity (rd,i6o) f rom Equation (l| 
for each line of sight. The observations reveal moderate 
FIR opacities toward Orion S and along the molecular ridge 
(Td,i 60 - 0.2). The FIR opacity then decreases from the Bar 
(Td,i 60 ^ 0.1) to the extended regions (r^o - 0.01). Note 
that the 160yum dus t emission toward Ori on BN/KL isopti- 
cally thick (Td,i 60 >l>[Goicoechea et al.[2015] and Sect. |5. 4.1). 

The total gas mass in the mapped region is computed from 
the optically thin dust emission at 250 pm using Equ ation ([2]). 
For the assumed grain properties (see Section |3.2| ), we de¬ 
rive M Gas , Total - 2600 Mq in the -0.9 pc x 1.4 pc mapped area 
(excluding Orion BN/KL that adds -100 M 0 ; |Goicoechea| 
et al.| 2015| ). This is equivalent to a surface density of 
^Gas-2000 Mq pc -2 . We note that a Rgd /^250 ratio higher than 
the assumed in this work will increase our gas mass estimate 
and the An-to-Ay conversion value accordingly. Unfortu¬ 
nately, the exact grain composition and optical properties in 


molecular clouds l ike OMC 1 are not fully constrained (Gold 
smith et al. T997| , and they likely change on small scales 
(Ara b et al.|2012| ). Owing to the similar gas masses re ported 
from large- scale^ 3 CO and submm continuu m mapping (Bally 


|et al.|1987l|Lis et al.|199 8 [ Bern e et al.|2014 ), we estimate that 
our M Gas and Ah values are uncertain by a factor of <2. 

3.4.2. L[C 11 7/LpiR and L[C 11 ]/L C o 2 - 1 Maps 

Figure [8jd) shows the spatial distribution of the L[C ii]/Lfir 
luminosity ratio toward OMC 1. L[C ii]/Lfir varies from 
the more extended and translucent cloud component, where 
Ah - 10 21-22 cm -2 and L[Cn] carries up to -1-5% of 
the FIR luminosity, to the dense molecular ridge, where 
A h - 10 24 cm -2 and L[Cn] only carries -0.01%. The mean 
(median) value in the map is L[Cii]/L F ir = 3.8 (2.9) x 10 -3 . 
This is similar to t he value observed toward translucent 
clouds of the Galaxy ( [Ingalls et al.|2002| ). Indeed, low angu¬ 
lar resolution surveys of the Milky Way show that the typ¬ 
ical L[Cii]/Lfir ratio in regions of weak [Cn] emission is 
~4xl0 -3 , and decreases to ~10 -3 to ward star-forming region s 
where the [Cn] line is brighter ( e.g.jNakagaw a et al.|1998| ). 

The Lfir map per pixel in Figure j8ja) can be used to Esti¬ 
mate, to first order, the spatial distribution of Go- Assuming 
that at large scales the entire column density of dust is heated 
by absorption of FUV photons, dust grains re-radiate in the 
FIR at a characteristic temperature (~ L^) and; 


~ 1 /fir( ergs ‘cm 2 sr ‘) 
0 ~ 2 1.3 x 10- 4 


(3) 


Line luminosity per pixel integrating all [C 11 ] spectral components. 


(see Hollenbach & Tiele ns[1999j ) where Gq is the FUV field 
in Habing unit^ The factor of 1/2 approximately takes into 

10 Using Ah (cm _2 /pixel) = A(H) + 2 A(H 2 ) = MG^/ipm^D 2 O), with 
Afpas obtained from Equation (2). 

11 At the distance of OMC 1 and considering the pixel size in Fig. [8(a), 
the Lfir scale per pixel is equivalent to Gq « 200 Lfir/Ld (per pixel). This 
relation neglects dust heating by embedded sources. Here we assume that 
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Figure 8. Luminosity and dust opacity per pixel maps toward OMC 1 convolved to 25" resolution, (a) Lfir map (Go ~ 200 Lfir/L© per pixel, see footnote 11), 
(b) L[Cn], (c) Leo 2-1 (from |Berne et al.|2014|> , (d) L[Cii]/Lfir in %, (e) L[Cii]/Lco 2 -i, and (f) dust opacity at 160 jum from SED fits. The white squares in 
panel (c) represen t pos itions of low Lco 2 -i and high L[Cii]/Lco 2 -i ratios (up to -5000) in panel (e) see also Figure fiTf. The different circles and regions are 
defined in Section [T2] The black squares indicate the positions of the spectra shown in Figure [3] The white square areas in the central BN/KL region indicate 
positions where the Herschel photometric observations are saturated. 


account the absorption of visible photons by dust. The max¬ 
imum value near the Trapezium is Go-10 5 , whereas the low¬ 
est values in the mapped area are Go >5 00. The mean value 
inside the 200" radius circle is Go-2xl0 4 (the Go used in 
the representative PDR models shown in Figure [5}. These 
FUV fluxes are in good agreement with early estimates based 
on lo wer angular resolution observations and PDR mo delling 


(e.g., |Tielens & Hollenbach| 1985b Stacey et al.|1993| ). 

The [Cn] and CO 2-1 luminosities in the observed area 
are L[Cn]=173L 0 and L c o2-i= 0.15L o . Figure [ 8 je) shows 
the L[Cii]/L c0 2 -i map. The spatial distribution is markedly 
different from that of L[Cii]/Lfir. In particular, the 
L[Cii]/L c0 2 -i ratio shows differences locally, from a few 
hundred in the strong outflows around Orion BN/KL, to 
-5000 in the peculiar northeastern positions almost devoid 
of CO emission (see white marks in the Leo 2-1 map shown in 
Figure [ 8 jc)). At other positions of the map, the ratio reaches 

owing to the high FUV field in the region, this extra heating does not domi¬ 
nate at large scales. For the same reason, and owing to the warm dust grain 
temperatures, we neglect the dust emission at >500 /urn in the Go estimation. 


high values (> 1000 ), either because the observations trace re¬ 
gions of low extinction (i.e., low CO columns) or nearly edge- 
on PDRs for which the C + /CO transition is spatially resolved. 
Table [2] summarizes the mean values extracted from the maps 
toward the different regions of the cloud. 

4. ANALYSIS 

4.1. [Cn] 158 fim Excitation and C + Column Densities 

Collisions with electrons, H atoms and H 2 molecules 
contribute to the excitation of the [Cn]^/^/^- 2 /^ fine 


coniriDuie to me excitation oi m e [1.11] r 3 / 2 - *1/2 lin e 
structure levels (see discussion by Goldsmith et al. 201 2| ) . 
The critical densities (n CY ) for collisions of C + with, e. 


1 ‘ — w V'wry ■ vymjiViio yi. _ ^ 

Hand H 2 are -44 e/cm 3 at 77 =8000 K (jWilson & Bell 


2002), and -3000 (2400) H/cm 3 ( [Barinovs e t al. 2005 )) and 
-4500 (3400) H 2 / cm 3 at T k = 100 (500) K ( Wiesenf eld & 
Goldsmit h|2014| ). Owing to the similar collisional rates with 
H 2 and H, the [C 11 ] excitation in neutral gas is nearly indepen¬ 
dent of the molecular fraction. As the gas density increases 
above n cv , the excitation temperature L ex of the 3 L ) 3 / 2 - 2 L > i / 2 
transition tends to thermalize with T k . For optically thick 
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Table 2 

Average physical parameters toward different regions of OMC 1 



L[C u]/L F ir 

£co2-iTfir 

L[Cii]/Lco2-i 

LpiR/Moas 

G b 

^d,160 





(Lq/Mq) 

(Habing field) 


OMC 1. 

3.8(5.4)xl0 -3 

3.6(6.2)xl0 -6 

1200(900) 

43(27) 

7.5X10 3 

0.03(0.03) 

(full map a ) 

Extended Cloud. 

4.7(6.0)xl0 -3 

4.5(7.5)xl0 -6 

1400(900) 

35(20) 

3.2xl0 3 

0.02(0.02) 

(F>0.4pc) 

Dense PDR interface. 

l.l(0.5)xlCr 3 

6.6(2.2)xl0 -7 

1600(500) 

87(24) 

3.3xl0 4 

0.08(0.05) 

(0.16<R<0.24pc annulus a ) 
Toward H 11 around Trapezium .. 
(R<0.1pc) 

9.3(3.9)xl0 -4 

5.4(1.8)xl0 -7 

1700(400) 

126(18) 

3.3xl0 4 

0.05(0.02) 


Note. — Average values over the different template cloud regions defined in Section |3.2| Numbers in parenthesis represent the standard deviation (lcr) of 
each sample. Lfir refers to the 40-500 pm range. The central Orion BN/KL region (50"x50") is not included in the computation of the average values due 
to saturation of the photometric measurements. A low L\C ii]/Lfir-(2-8)x 10 -5 luminosity ratio and 160-1-3 was inferred from HerschelfPACS spectroscopic 
observations jGoicoechea et al .|2015) . The mean L[C ii]/Lco2-i ratio in this region is ~270. ^Estimated from /fir and Equation]^] 


lines, line-trapping further reduces the critical densities, thus 
reaching thermalization (r ex =7k) at lower densities. Also for 
optically thick lines, a strong FIR background with high Tbg 
values increases T ex but reduces the [C n] line intensity. 

When nu ^ the [C n] excitation becomes sub-thermal 
(T ex < 7k). If significant C + columns of low-density gas exist 
in front of a FIR continuum source, the [C ii] line will be seen 
in absorption if T ex < k(1+ 9 9 1 1 2 2/r - (e.g.,|Gerin et al. 20T5| , 
where hv/k = 91.2 K is the equivalent temperature at 


1900.537 GHz, and T c = J(T hg ) = 


91.2 


is the continuum 


e 9L2/T ^-l 

brightness temperature at 158 /mi. [Cn] line self-absorptions 
can also be produced by excitation gradients within the [C n] 


emitting region (e.g., |Plume et al.|2004||Graf et al.|2012| ). 
One potential problem in interpreting the [Cn] 158 yt/m 


emission is the difficulty inferring the line opacity. Toward 
many positions, our observations resolve the [ 13 Cn] F= 2-1 
linqT^J arising from the main spectral component at 
vlsr - +9.5 km s -1 (see spectra in Figure 3] and maps in the 
highest-velocity panels of Figures [6] anam). The second 
brightest component, the [ 13 Cn] F=L0 line, is only detected 
toward a few positions of very bright [C n] emission (e.g. the 
Orion Bar and BN/KL, see Figure [5]). Assuming that the to¬ 
tal integrated intensity of the threeT 3 Cn] lines is 7[ 13 Cii] = 
h= 2-i/0.625, one can derive the [Cn] line opacity froirpj 


1 - 0.625 r P ([ 12 CII])/r P ([ i3 CII], F = 2-1) 


T [CII] 


[ 12 C/ 13 C] 


(4) 


(e.g.,|Ossenkopf et al. 2013J) where T[cii] is the [Cn] 158 /im 
opacity at line center, Y’ P is the line peak main beam tem¬ 
perature (continuum-su btracted) and [ 12 C/ 13 C ] is the isotopic 
ratio (~67 in Orion, Langer & Penzias||1990 ). For optically 
thick [Cn] 158yum emission, the excitation temperature (T ex ) 


12 Corrected line strengths for the [ 13 C n] F=2-l, 1-0 and 1-1 hyperfine 
lines at 1900.466, 1900.950 and 1900.136 GHz respectively have been pub¬ 
lished by Ossenkopf et al. ( 2013!. The relative line strengt hs are 0 . 625, 0 .250 
and 0.125. these are slightly differe nt from those of |Cooksy et a l. (1986) that 
were used in the analysis of jStacey eta l. V 199Tb) and fBoreiKO & Betz (19 96), 
and resulted in 12 C/ 13 C isotopic ratios lower than the usual value of ~67. 

13 Toward the Orion Bar and BN/KL, the [ 13 Cn] F=2-l line is blended 
with the red-shifted [C ii] emission. To avoid overestimating the [C n] 158 pm 
line opacity toward these positions, we used the [ 13 C ii] F=l-0 line (which is 
not blended with any spectral feature) and assumed /[ 13 Cn] = /p =1 _o/0.250. 


of the 2 P 3/2 - 2 Pi /2 transition can be computed from: 

91.2 

= 771-513-7 (K)- ( 5 ) 

ln(l + rp( i2 c+>+/(r bg )) 

At these FIR wavelengths, the background is dominated by 
dust emission so that the continuum level adjacent to the [C n] 
line, for an equivalent background temperature of 20-35 K, is 
T c = 1 - 7 K, roughly the range of continuum levels detected 
by HIFI. For the main [Cn] 158 yum spectral component we 
find r P ( 12 C + ) » T c , thus the continuum emission has a minor 
effect on T QX . However, the situation can be different in more 
extreme and luminous environments. 

Once T ex ( 2 P 3/2 - 2 P 1 / 2 ) and T[cn] have been determined, the 
C + column density can be computed from: 

1 + 2 e ~ 9L2/T ™ 

N(C + ) - 1.5x 10 17 - _ ^_ 912IT ^ t [cii] • Avfwhm (cm -2 ) 

( 6 ) 

where Avfwhm is the line width in kms 1 of an as¬ 
sumed Gaussian line-profile. From the [Cii] and [ 13 Cn] 
spectra averaged over the mapped region (main spectral 
component) we find < T ex 100 K, < T[cn] >-1.3 and 
< N(C + ) >^3x 10 18 cm -2 . The highest C + column den¬ 
sity peaks are found near the Trapezium, the Orion Bar, the 
East PDR and northeast of Orion BN/KL. In these emis¬ 
sion peaks, the excitation temperatures and line opacities 
(T ex = 250 - 300 K and T[cii] -1-2) are higher than the av¬ 
erage. The 25" beam-averaged C + column density toward 
these peaks is N(C + ) - 10 19 cm -2 . 

The inferred excitation temperatures can be used to fur¬ 
ther constrain the origin of the [Cn] emission at the 
FUV-illuminated surface of the molecular cloud. In par¬ 
ticular, the mean T ex >100K values inferred from the 
map require the fractional population of the 2 P 3 /2 level, 
ni 2 P 3 / 2 )/n( 2 P 1 / 2 + 2 P 3 / 2 ), to be >0.45. The representative PDR 
model in Figure[5]shows that the conditions for producing sig¬ 
nificant [C 11 ] emission are quite restrictive. For high FUV ra¬ 
diation fields, the C + layers extend to A v -3-4 mag but the 
gas temperature falls much faster for any reasonable den¬ 
sity. Figure [5] shows that a fractional population >0.45 of the 
2 P 3 /2 level requires high gas densities (^h-5x 10 3 for even 
r k ^500K, see Figure (10|). Hence, the bright [Cn] emission 
from OMC 1 traces narrow cloud surface layers (Ay<4 mag) 
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Figure 9. Column density per pixel maps toward Orion, (a) 7V(C + ) blue extracted from the blue-shifted emission component from vlsr= - 2 to +5 km -1 , 
(b) iV(C + ) 0MC1 extracted from the main emission component from vlsr-+5 to +17 km -1 , (c) N(C + ) Total =N(C + ) OMCl +N(C + ) hlue +N(C + ) < ~ 2km ^, where 
A/'(C + ) <-2km/s =10 17 cm -2 is adopted in all positions (see text), (d) N^ tal =N(R)+2N(R2) from the dust emission, (e) [A^(C + ) Total / J /VH° tal ]/1.4xlO -4 color map, 
this is roughly the fraction (/) of the total column density traced by [C n] toward Orion. This fraction scales with the total column density of dust a long each line 
of sight. The 160gm dust opacity iso-contours at Td,i60=0.1 (inner) and 0.05 (outer) are shown in black.The different circles are defined in Section|T2] 


and not the entire line of sight. 

In the extended regions, the [ 13 Cn] lines are not detected 
and the [Cn] 158 /im line is weaker. This su ggests that the 
line is optically thin or effectively thin (Goldsmith et al. 


|2012| ). The line opacity can be derived if one knows T ex . 
Here we make the approximation T[cn] - 7p([C n])/J(T QX ). 
Figure [9ja) shows a column density map, computed under 
these assumptions, of the blue-shifted component (vlsr- - 2 
to +5 kms -1 ) that includes the Dark Bay and the North¬ 
ern Dark Lane emission. We adopted a constant T ex ^80K 
(intermediate between 7\ in the OMCl/Hn interfaces and 
that in the atomic Veil ( [Abel et aLl|2009| ). The resulting 
N(C + ) toward the Dark Bay (^1.4xlO i0 cm -2 ) is equivalent 
to Ay^/V(C + )/4xl0 17 ^3.5 mag (i.e., significant H 2 must be 
present) assuming that the bulk of the gas-phase carbon is in 
C + , with an abundance of 1.4xl0 -4 ( Sofia et al.|2004 ). This 
value is consistent with the extinction peak determined toward 
the Dark Bay ( [O’Dell & Yusef-Zadeh 2000| ). Other sightlines 
with V(C + )<8xl0 17 cm -2 (Ay<2mag) are likely associated 
with C + in nearly atomic H i gas (e.g., toward the Traprezium). 
The integrated C + column density in the map (vlsr--2 to 
+5 kms -1 component) is compatible w ith t he observed [Cn] 
fractional luminosity of ~12% (Section [33] ). 

In addition, Figure |9jb) shows a N(C + ) map of the OMC 1 
component (vlsr-5 to +17 km s -1 ) adopting r ex ^100 K from 


the map averaged spectrum. Figure [9fc) shows the resulting 
map of the total C + column density along each line of sight, 
V(C + ) Total =V(C + ) 0MC1 +N(C + ) blue +V(C + ) <-2 , where we 
have added, at each position, A(C + ) <-2km/s =10 17 cm -2 
from additional foreground emission components at 
vlsr< - 2 kms -1 . The integrated column density in the 
<-2 kms -1 component is consistent with the ~3% contri¬ 
bution to the total [Cn] luminosity. This is equivalent to 
A<- 2 km/s^0 25 mag (i.e., fully atomic or ionized gas). We 
finally note that along many sightlines, the C + column in the 
Veil is comparable to that computed in the OMC 1 surface. In 
these approximate maps of the C + components without CO 
counterpart, ~30 % of the C + mass is in CO-dark H 2 gas. 

Figure |9^d) shows the V^ otal column density map ob¬ 
tained from the dust emission 10 . Both A(C + ) Total and 
Vjjj 0tal maps can be used to estimate the fraction (/) of 
the total column density traced by [Cn] along each line of 
sight. In particular, Figure [9je) shows the distribution of 
the /=[V(C + )/Vh]/1.4x10 -4 ratio at 25" resolution. Note 
that this map is similar to the L[Cii]/Lfir map shown in 
Figure [8jd). The lowest values / <$c 1 are inferred toward the 
high Nu column density ridge (bluish regions with / <5 %). 
These fractions are consistent with the presence of a massive 
molecular cloud behind the [C n] emitting surfaces. The typ- 
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^ocai We used a constant col- 
the mean value we in- 


ical fractions outside the ridge are / ^20-30 % (greenish re¬ 
gions). Far from the Trapezium, toward the extended cloud 
component, fractions higher than />50 % are only inferred 
toward the most translucent lines of sight. The [C n] emis¬ 
sion toward the dense PDRs around the Trapezium is slightly 
optically thick and T ex is higher than the average value as¬ 
sumed in the / map. The corrected fraction using the column 
of V(C + )^10 19 cm -2 computed from [ 13 C n] is 40 %. 

We also used the C + column map to compute the “PDR 
mass” in OMC1. In particular, we integrate the equiva¬ 
lent Vh,pdr=^(C + ) omc1 /1-4x10 -4 column traced by [Cn] as 
M PDR - yumn Z* C^hpdr^')’ where A t is the area of pixel i, 
mu is the H atom mass and p is the mean atomic weight. We 
derive M PDR ^200M o , i.e., the PDR mass fraction traced by 
[C n] is -8 % (within a factor of -2). 

4.1.1. [Cn] Radiative Transfer Models 

In order to constrain the range of physical conditions 
that reproduce the observed [Cn] 158/mi emission to¬ 
ward OMC 1, we have run a grid of non-local and non- 
LTE radiative transfer models of the 2 P ^/2 ~ 2 P 1/2 doublet 
that take into account line trapping and opa city broaden¬ 
ing (see Appendix in |Goicoechea et al.|[2006| ). For consis¬ 
tency, we included illumination by a FIR background contin¬ 
uum field and a combined treatment of gas and dust emis¬ 
sion/absorption. Hence, at a given wavelength the opac- 
ity is given by t a = r [C n] + r^. 

umn density of N( C + )=3xl0^ cm" 
fer from the OMC 1 map, and a nonthermal velocity dis¬ 
persion cr nt h = 1.3 kms" 1 (see Sec. 3.3.2 ). The models use 
the latest available collisional rates of 2 + with e , H and 
H 2 (with a LTE H 2 ortho-to-para ratio at each 7k) and as¬ 
sume n e = [C + /H] n H = 1.4 x 10 -4 n H and n(H)=0.25 n(H 2 ) 
(a molecular fraction, 2 n(H 2 )/[n(il)+2 n(H 2 )], of ^0.9). For 
gas densities ^h=2x10 5 cm -3 and 7k=500 K, this choice of 
parameters implies that the [C 11 ] 2 P3/2~ 2 Pl/2 collisional rate is 
due to collisions with H 2 molecules (-70%), H atoms (-25%) 
and electrons (-5%). The FIR background radiation field is 
modelled as thermal emission at 30 K (roughly the mean dust 
temperature we derive from SED fits) and the temperature of 
the dust grains mixed with the C + atoms is set to 50 K (typical 
of highly illuminated PDRs at A v « 2 mag). 

Figure [lO] shows model results in the form of iso-T ex 
contours. The mean [Cn] excitation temperature in the 
Herschel/ HIFI map is T ex ^100 K. This plot shows that the 
minimum gas density required to produce this T ex =100K 
is ^h-5x10 3 cm -3 . Toward the Hn/OMC 1 interfaces (the 
spherical shell of PDRs around the Trapezium) T ex reaches 
higher values (-250-300 K). This sets a higher limit to the gas 
density of tth^IO 5 cm -3 (including the Orion Bar). The more 
elevated temperatures in these dense PDRs (7k>300K) are 
consistent with the 7k=400-700K values inferred from H 2 
pure rotational lines ( [Allers et aL i; 2005| and with T k ^540 K 
derived from H i emission, both measur ements toward the 
Orion Bar PDR |van der Werf et al.l2013[ ). 

Finally, the presence of a foreground halo of low-excitation 
material can produce [C 11 ] absorption. Such a halo will also 
absorb or scatter the [Cn] line photons emitted inside the 
cloud, if the absorption is produced at similar velocities. This 
may be the case of components A and B of the Hi absorp- 
tion in the Veil. These components peak at vlsr- 0-6 kms -1 
( [van der Werf et al.|2013| ) and coincide with the blue-wing 
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Figure 10. Synthetic excitation temperatures for the [C n] 158 fim transition 
obtained from a grid of nonlocal non-LTE models (see text). 

of the [Cn] emission seen in the main spectral component. 
The gas temperature and density in components A and B of 
Orion’s Veil are estimated from Hi lines to 77-50-80 K and 
«(H)^10 2 5-3 ' 4 cm" 3 dTroland et al.|l989l|Abel et al.|2006b| l. 
For the FIR background in OMC 1 and for the parameters 
in our grid, we do predict [Cn] absorption for low excita¬ 
tion gas (white line in Figure [T0|. Hence, the [Cn] emis¬ 
sion dips seen in some velocity intervals may be produced by 
[Cn] absorption in the Veil (e.g. toward the Orion Bar PDR 
at vlsr-+1 km s -1 , toward the southwest (-150",-300") posi¬ 
tion at Vlsr-+3 km s -1 or toward the (-62",250") position in 
the northern streamer at VLSR-+6km s _1 , see channel maps). 

In any case, our [C n] excitation models confirm that the 
gas density in the FUV-illuminated face of OMC 1 is high. 
The blue-shifted [C n] emission caused by the Veil also shows 
that the density must be relatively high in these foreground 
components. This agrees with the general notion that the el¬ 
evated stellar density in the Orion cluster is a consequence of 
hundreds of st ars that formed from a remarkably dense co n- 
densation (e.g., [Genzel & Stutzki|1989t|Rivilla et al.|2013| ). 


4.2. Correlation Diagrams 

Figure |TT] shows correlation plots between several quanti¬ 
ties that are also used in the discussion of the [Cn] extra- 
galactic emission (one point, one galaxy). Here we present 
their spatial correlation in a template massive star-forming re¬ 
gion. In these plots, the black points represent positions out¬ 
side the 200" (0.4 pc) radius circle, the “extended cloud face” 
component. Green squares represent lines of sight toward the 
“Hu region” component around the Trapezium stars, and red 
triangles represent positions in the “dense PDR interface” an¬ 
nulus. Gray points are other positions inside the 200" circle 
that do not fall in the latter two categories. 

Figure [Tlja) shows that the [Cn]/FIR ratio decreases 
with increasing FIR intensity. In OMC 1 we find that 
log ([C n]/FIR)^-0.6-log (7piR[erg s -1 cm -2 sr _1 ])-2.6, with a 
correlation coefficient of p=-0.93. This is reminiscent of 
the “[C n] deficit” observed in local ULIRGs when compared 



curs at orders of magnitude lower FIR luminosities and seems 
dominated by the strong variations of 7 fir in the region. The 
blue line in Fig. [lTJa) shows a linear model for a constant 
/ cloua ( 


[C n] intensity < 


[ (the median value of the map sample). 
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Figure 11. Correlation plots extracted from the OMC 1 maps and convolved to a uniform 25" resolution. Black points represent positions outside the 
200" (0.4pc) radius circle (“extended cloud face” component). Green squares represent positions toward the Hu region component around the Trapezium. 
Red triangles represent positions in the “dense PDR interface” annulus (see e.g., Figurep). Gray points are other positions inside the 200" (0.4 pc) radius circle 
that do not fall in the latter two categories. Two models (blue lines) are plotted in panels (a) and (b), one for a constant [Cn] intensity cloud (a) and another for a 
cloud of constant T& (b). Both models are defined to intercept the median L[Cii]/Lfir, /fir and 7^,160 of the map (cyan dots). 

tween dust opacity and FIR intensity shown in Figure JTTJb). 
However, the inferred dust temperatures toward some lines of 
sight of small Td,i 60 are lower than 30 K, resulting in lower 
Lfir values than those expected in a constant Tj cloud (blue 
line Figure pTJb)), as well as in high L [C ii]/Lfir ratios in Fig¬ 
ure [Tlja). On the other hand, the dust is warmer toward the 
dense PDRs and H n region, and high / FIR intensities are ob¬ 
served above the constant- Tj blue line in Figure p~Tl(b). The 
higher dust temperatures toward the central of OMC 1 likely 
reflect higher Go values and also contribution from internal 
dust heating by embedded sources. However, their associ¬ 
ated L[C ii]/L fir ratios lie above the linear model that assumes 
constant [C n] emission in Figure |TTJa), revealing enhanced 
[C n] intensity toward these regions. Still, the FIR intensity 
increases more than [Cn] (Figure pT^d)), and the resulting 
L[C ii]/L fir ratios are low toward lines of sight of high dust 


Figu re pT^ b) shows that / F i R increases with Td,i 60 - 
Figure pT^c) shows a L[Cn]/L F i R -7d scatter plot. In many 
studies, in which a complete SED cannot be constructed, T d 
is replaced by the F^o/Fioo or Fjo/F^q flux ratio. Theo¬ 
retical models show that the photoelectic heating efficiency 
decreases for high GoM h values (also producing higher dust 
temperatures), because grains and PAHs get more positively 
charged and it beco mes more difficult to photo- eject addi¬ 
tional electrons ( Tielens & Hollenbach| 1985) |Okada et al.| 
|2013| ). As the grain charge increases for high Go fields, the gas 
heating and cooling are reduced and a larger fraction of the 
FUV field is converted into dust FIR luminosity (thus lower 
L[Cii]/L fir ratios can be expected). Although such a trend 
is suggested by Figure JTTJc), the dust temperatures obtained 
from the SED fits, 7 o= 3T±4 K (lcr), do not change much over 
the mapped region and other effects seem to dominate the ob¬ 
served L[C ii]/L F i R variations at large spatial scales. 

In particular, most of the dust emission in the extended 
cloud is consistent with an roughly constant temperature of 
Td-30 K. This leads to the approximately linear relation be- 


temperatures, see Figure [TTJc) and Section [5Al 
Figure |TTJd) shows that, except for high FIR in¬ 
tensities, the [Cn]158yum intensities are reasonably 
well correlated with the FIR emission (p=0.75), with 
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I[ cii]-4.1x10 _4 -/fir+2.2x 10“ 3 (in ergs -1 cm -2 sr -1 units). 
The [Cn] intensities, however, are better correlated with 
the PAH emission traced by the IRAC8/im image (p= 0.91, 
Figure [lTJf)), with / [C ii]-2.6xl0 _2 -/iR A c8 y um+1.6xl0“ 3 (also 
in erg s -1 cm -2 sr -1 units). As opposed to FIR, this correlation 
seems to hold for high [Cn] and PAH intensities. Finally, 
Figure [Tlje) shows a correlation plot of L [C ii]/Lfir with 
LpiR/Mcas- This plot shows much more scatter, but suggests 
that L[Cii]/L fir decreases from the extended cloud to the 
“Hu region” component (highest LpiR/Moas values). 

5. DISCUSSION 

5.1. [Cn], CO and FIR Diagnostic Power 

Dense PDR models predict that the low -J 12 CO lines sat¬ 
urate (r»l) close to the CO abundance peak and thus the 
L[C n]/Lco ratio depends mostly on N(C + ) and 7\, i.e., on Go 
and nu (see e.g., PDR models by [[Kaufman et al.|1999| . For 
densities below ft cr ~ 10 4 cm -3 an d Gp<10 3 , the [Cn] 15 8 /am 
line dominates the gas cooling (Kaufm an et al.|T999| ) and 
the L[C ii]/L fir ratio provides a good approximation"to the 
grain photoelectric heating efficiency. For higher densities or 
Go values (i.e., the OMC 1/Hn interfaces around the Trapez¬ 
ium), [Cn] collisional de-excitation plays a role and [Oi] 
starts to dominate the warm gas cooling. [Oi] 63,145yum 
fines are certainly useful ( Herrmann et al. 1991) although they 
are also b right in dissociative sho cks where the [C n] fine is 
faint (e.g., Goicoechea et al. |2015| ). In addition, the opacity of 
the [Oi] 63//m line is higher (t[oi] = 1.6t[cii] for a [0/C]=2.2 
abundance ratio) and lies at wavelengths of higher dust opac¬ 
ity (2.5-6.3 times higher for /3=l-2). The [Oi]63yum fine 
opacity problem is difficult to circumvent, because 18 0 does 
not have nuclear spin (producing hyperfine splittings in 13 C + ), 
and the isotopic shift with respect to 16 0 is only ~1.5 kms -1 
(+23.5MHz, Brown et al.|1993 ). Therefore, [Cn] is a robust 
tracer of the presence of FUV photons also in dense gas. 

As noted in the literature, L[Cii]/L fir approximately fol¬ 
lows L C o 2 -i/T F ir (Figure [12]). For most of the mapped 
positions, the observed luminosi ties appear we ll correlated 
(p= 0.89 and slope of ^775). [Stacey et al. \\99l) sug¬ 
gested that, to a first approximation, the L[Cii]/L fir and 
Tco/T F ir luminosity ratios can be compared with PDR mod¬ 
els to constrain the FUV radiation field and the gas density. 
Most of the points in Figure [12] roughly fie on a straight 
line. In the frame of PDR models (e.g., Kaufman et al. 
1999) this line approximately follows the decrease of Go 
(from left to right) and the decrease of nu (from bottom 
to top). Comparing our spatially resolved l uminosity ra- 


10 1 - 


tio|^| (Figure [l5]) with the plot in Figure 5 of 


_ : Stacey et al. 

(2010), based on PDR models of Kaufm arTet aL| (\99Vf for a 
Ay—10 mag slab, one sees that the lowest L[C ii]/L F ir^3x 10 -4 
and Lco2-i/T F ir-3x 10 -6 values near the Trapezium are con¬ 
sistent with high Go>10 4 fields and tth^IO 5 cm -3 . On the 
other hand, the highest L[Cii]/L F ir> 10 -2 values, correspond¬ 
ing to Lco 2 -i/T F ir< 10 -4 ratios in the extended cloud compo¬ 
nent are consistent with lower FUV fields and gas densities, 
Go>10 2 and ^h- 10 3-4 cm -3 respectively. 

In addition, our maps reveal regions in which the luminosity 
ratios lie outside the parameter space allowed by Ay =10 mag 
PDR models (blue pentagons in Figure [12]). These are spe¬ 
cific positions with very high L[Cii]/Lco2-i- 5000 luminosi- 

14 Scaling our data by Lfir(40-5005 ^ 42 - 122 /im and assuming 
4:02-1-8 Ico 1-0 (erg s -1 cm -2 sr -1 ) for optically thick thermalized CO gas. 


J 1 

u 


10 r 


1—1—1—1 1 1 11 1 -1—1—1—1 1 111 1 -1—1—1—1 1 111 1 

Extended 0MC1 face ° 

Interface, dense PDRs a 0 ° 0 o ° 

HII around Trapezium u v ° 8 °°°° 
Low L co NE—region 


10 r 



10 ° 10 
Lco 2 -i/LfIR 

Figure 12. Same as Fig. JIT] but for L[Cii]/Lfir versus Lco2-i/Lfir- The 
blue pentagons represent lines of sight of low L(CO) and high L[C ii]/Lfir 
ratios (up to -5000) at the northeast of the map, see Figures Me) and[8je). 
Arrows show the approximate dire ctions of increasing Gp and increasing n u 
expected from PDR models (e.g., | Kaufman et al.|1999[|Stacey et al.|2010) . 

ties and low Td,i 60 opacities at the northeast of the map. Since 
they do not coincide with the Hu region, we conclude that 
they must be associated with lower density neutral gas at low 
Ay (thus low CO abundances for Go»l values). Indeed, a 
much better match with PDR models is found by considering 
the Ay ^3 mag cloud models of Kauf man et al. ( 1999) . 


5.2. [Cn] and the Star Formation Rate in Orion 

Herschel/YHFI has allowed the first velocity-resolved [C 11 ] 
survey of the Galactic plane, in which the contribution of 
the different ISM pha ses t o the [C 11 ] emission has been con¬ 
strained (e.g., [Pineda et al.|2013| . The combined [Cn] emis¬ 
sion from all components added together correlates with the 
SFR in similar manner as in nearby galaxies. Including obser¬ 
vations of the Milky Way, the LMC and several n earby g alax- 
ies dominated by star formation, |Pineda et al.| (|2014j) find: 
log (SFR [M 0 yr -1 ])^0.89 log (L[C ii][er g r i ])-36.3. Using 
this correlation and the total [C 11 ] luminosity toward 
OMC1, we derive SFR^3.8xl0 -5 M 0 yr -1 in the mapped 
area. This is very similar to the rate_ obtained using 
SFR(M 0 yr -1 A Lfir/5.8 x 10 9 L o (e.g., |Kennicutt| l998) and 
the observed luminosity toward OMC 1 (Lfir-2x 10 5 L q ). 

In addition, extinction maps of galactic star-forming 
regions suggest that the SFR correlates with the 
mass (M r ) of the UV-shielded (Ay>7 mag) dense gas 
(ft(H 2 )> 10 4 cm~ 3 ), the reservoir that ultimatel y forms stars 
(e.g., Andre et al.|2010[|Lada et al.|2010| | 20l2| . These 
studies suggest SFR [M e yr -1 ] ^ (1.2-1.8)xl0 -8 M' (M 0 ), 
and this corrrelation holds from star-forming r egions to 
galaxies ( [Gao & Solomon||2004[ |Wu et al.]|2005| ). Taking 
the masses inferred from our observations, we estimate 
M'— AfGas,Total pdr— 2400 M© (using the total mass deduced 
from the dust emission and subtracting the PDR contribution 
inferred from [Cn]), and use the above relation to derive 
SFR^(2.9-4.3)xlO -5 M 0 yr _1 . Therefore, the different 
methods provide consistent SFR values in Orion. 

Taking into account the mapped area, our derived 
SFR is equivalent to a very high SFR surface density 
of S SF r=(2.3-3.4)x 10 -5 M 0 yr _1 pc -2 . Together with the 
high surface density inferred from the mass estimates 
(^Gas-2000 M q pc -2 ), these values place the core of Orion at 
an extreme of a Kennicutt-Schmidt-type surface density rela¬ 
tion (2 s F r oc with a close to 1). These, and even higher 
surface densities are also inferred toward lumin ous g alaxy 
mergers hosting intense star formation ( [Genzel et al.|2010| ). 
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Figure 13. Spatial distribution of the 1 -fir/^G as ratio in units of L© M~ x 
per pixel (colored map) and H41a integrated intensity map in gray contours 
(from 2 to 22 Kkms -1 in steps of 4 Kkms -1 ). The white thick curve shows 
the l-FiR/MGas=80 Lq Mq 1 iso-contour. Note the good correspondence be¬ 
tween the ionized gas emission and the regions of high 2-fir /Mjas ratios. 
The red, yellow and the gray stars show the position of the Orion hot core, 
Orion S and the Trapezium stars respectively. The black squares show the 
positions of the spectra shown in Figure[3] 

5.3. L[C ii7 /Z-fir Variations and Extragalactic Link 

The reduced L[C ii]/Lfir luminosity observed toward local 
ULIRGs (L F ir>10 12 Lq) has been difficult to interpret. Dust 
extinction, optically thick [C n] emission, reduced photoelec¬ 
tric heating efficiency or soft UV field s from less massive stars 
have been invoked (e.g., |Luhman et al. 1998| ). In addition, ex¬ 
tragalactic observations with Herschel suggest that the “line 
deficit” (relative to Lfir or Lpm/Moas) also applies to other 
FIR fine structure lines such as [Qi], [Nn] or [Om] arisin g 
either from PDRs or H n regions ( Gracia-Carpio et al.[201 1 ). 

The LpiR/Moas ratio in particular, is expected to be propor- 
tional to the ionization parameter U (the number of H ionizing 
photons, 2(H), divided by the density n h and by the distance 
to the ionizing stars). The FIR line deficit becomes more ap- 
paren t in galaxies with very high TpiR/Mcas^O Mp 1 val¬ 


ues ( [Gracia-Carpio et al.|201 l[|Gonzalez- Alfonso et al. 2015). 
This threshold has been suggested to separate normal galaxies 

from ultra-luminous mergers characterized by very high SFRs 
( [Genzel et al.|2010) . 

These authors suggest that the ISM conditions and the 
star-formation efficiency in these galaxies are different from 
those in normal and starburst galaxies. In this con¬ 
text, the FIR line deficits (including [Cii]) would be ex¬ 
plained by the much higher ionization parameters and 
more extreme conditions of th e Hn/PDR regions they host 
(e.g.,[Gracia-Carpio et al. 201 1| ). Photoionization models in¬ 
cluding dust grains indeed predict reduce d L[C ii] /L F ir values 

One of the 


for high log U > -2 parameters { [Abel et al.l! 
predicted effects is that grains absorb a higher fraction of the 
available UV photons relative to the gas, leaving fewer pho¬ 
tons to ionize and heat the gas. Models of the radiative pres¬ 
sure on grains do show that dusty H n regions have a different 


ionization balance and that they compress the gas and dust 
into a dense ionized shell ( |Draine|2011[ ). This is very sugges¬ 
tive of the shell-like structure revealed by the [C n] and H41cr 
emission around the Trapezium. On the other hand, |Draine 
P011| ) concludes that owing to the accumulation of gas and 
dust in a dense shell, the fraction of photons that ionize the 
gas actually increases (relative to a uniform Hu region). In 
this context, the effects of increasing U may not be so severe. 

Observational evidence for t he presence of PAH and grain s 
in Hu regions is increasing ( [Ochsendorf & Tielens|2015] ). 
Their role in modulating the FUV field reaching the neutral 
cloud surface needs to be taken into account. The H n region 
around the Trapezium contains dust and is indeed character¬ 
ized by a relatively high ionization parameter (log U- -1.48 
and <2(H)~10 49 s _1 , Baldwin et al.][1991). The presence of 
grains in the ionized nebula can be inferred from the depleted 
abundances of refractory elements_like Mg, Si and Fe from 
the gas phase (e.g., |Sim6n-Diaz & Stasinska|2011| . 

In Orion, the sightlines with LFiR/M Ga s>80 L© M~ l show a 
very good spatial correspondence with the dense H n ionized 
gas traced by the H41 a line and so they do trace regions of 
high ionization parameter (white thick contour in Figure [13]) 
In particular, the highest T F iR/M Ga s ratios are observed close 
to the H41 a emission peaks: the dense ionized gas toward the 
Trapezium cluster and also the ionization front that precedes 
the Orion Bar PDR. These regions are characterized by more 
extreme UV irradiation conditions (Go>10 4 ). 

We finally note that the luminosity of the large-scale vi- 
brationally excited H 2 line emission toward OMC 1 is domi¬ 
nated by FUV-exited fluorescent emission in PDR gas (>98%, 
see |Luhman et al.|l994 ) even though the collisionally excited 
(shocked) H 2 dominates the small-scale emission of high sur¬ 
face brightness (e.g., toward BN/KL and other outflows). Like 
the extended [C 11 ] emission we have mapped, this FUV-exited 
H 2 emission traces the FUV-illuminated face of Orion and 
may also dominate the extragalactic near-IR H 2 emission. 


Figure 14 


5.4. An Explanation for the 
shows 


“[Cn] Deficit” 

the spatial distribution of the two 
L[C nl/LpiR^egimes found toward OMC1. Cyan contours 
for low ~10 -4 -10 -3 ratios toward the core of the cloud, and 
red contours for high ~10 _3 -10 -2 ratios toward the extended 
cloud face component (increasing from the center of OMC 1 
and outwards). In Figure [l4| (left), the L[Cii]/Lfir contours 
are shown on top of the L^ir/Mq^ map. Figure [14] (right) 
shows the same L[C ii]/Lfir contours on top of the Td,i60 map. 

5.4.1. Geometry and N(C + ) Relative to Total Nu Column 

Since most of the [Cn] emission arises from Ay<4mag 
FUV-irradiated cloud surface layers, the local opacity of the 
dust grains mixed with the C + atoms (t^q) is generally 

small, - N h • pm^ Km/R g d ~ 0.001. Therefore, the 
dust opacity does not seem important locally, but it may be 
relevant if the [C 11 ] sources are embedded in large columns of 
dust or if many [C 11 ] emitting surfaces in a clumpy medium 
are included in the telescope beam. 

The average 160//m dust opacity inside the 200" circle re¬ 
gion surrounding the Trapezium is only moderate (r<i,i60 - 
0.06 or Ay>50mag). Therefore, the 160 pm continuum emis¬ 
sion traces the column density of dust in the molecular cloud 
that lies behind the ionized nebula, weighted by 7^. The [C 11 ] 
emission, however, does not arise from the entire line of sight 
but from the FUV-illuminated face of the cloud. 
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Figure 14. Spatial distribution of two regimes found for the L[Cii]/Lfir luminosity ratio. Cyan contours for low L[Cii]/Lfir ratios, 10 4 —1.5-10 3 toward the 
core of OMC 1, and red contours for high ratios, 2-10 _3 -10 -2 toward the extended cloud component (the L[C ii]/Lfir ratio increases from the center of the map 
and outwards), (a) The colored map shows the Lpm/Moas ratio, which to a first-approximation traces ionized gas and the ionization parameter U. (b) The colored 
map shows the spatial distribution of the dust opacity at 160 //m. The red, blue and the gray stars show the position of the Orion hot core, Orion S and the 
Trapezium stars respectively. The black squares denote the positions of the spectra shown in Figure[3] 


Figure p~5] shows the L[C ii]/L fir luminosity ratio as a func¬ 
tion of Td,i60- Both quantities are clearly tightly related. 
Therefore, the relative geometry of the illuminating stars 
with respect to the molecular cloud is important to under¬ 
stand the variations of L[Cii]/Lfir. Toward OMC 1, it can 
be idealized as a slab of constant [C n] emission arising from 
the foreground face of the molecular cloud (see a simpli¬ 
fied sketch in Figure [T6|z). In this case, and for constant r d , 
L[C ii]/Lfir would vary as 1/(1 -^ _Td160 ), with Td,i 60 being pro¬ 
portional to the dust column toward each line of sighp] When 



Figure 15. Observed L[C ii]/7,fir ratios as a function of the 160 //m dust 
opacity. The different template regions are defined in Section |T2| The in¬ 
set panel shows an extension to higher dust opacity positions toward BN/KL. 
The blue continuos, magenta dashed, and black dotted curves show models 
of constant 7[C n]/B(T d ) emission for different assumptions of the [C n] emis¬ 
sion location with respect to the cloud dust emission: [C n] in the foreground 
surface (as in OMC 1), mixed C + gas and grains, and embedded [C n] sources 
behind the cloud respectively (see sketch in Figure [76) . Models are defined 
to intercept the median L[C ii]/Afir and Td,i 60 valuesof the map (cyan dot). 


this slab model is defined to intercept the median L[C ii]/L F ir 
and Td,i 60 values of the sample, the resulting curve satisfacto¬ 
rily reproduces the ~3 orders of magnitude L[C ii]/Lfir varia¬ 
tions observed throughout the map (blue curve in Figure p~5). 
The observed scatter with respect to the model is likely pro¬ 
duced by small variations of T d , grain charge, and of I[C n] in 
lin es of si ght with the same r d j 60 but different Go. 

Gerin et ak| ( |2015| ) have carried out small maps of the 


[Cn] 158yum line toward a limited samplq ib | of more dis¬ 
tant massive star-forming regions in the Milky Way (W31C, 
W49N, W51, DR21(OH), ...). These observations show that 
the [Cn] absorption produced by diffuse clouds on the line 
of sight can reduce the L[Cii]/Lfir ratio by ~50 %, if the 
Vlsr range of the absorbing clouds is not spectrally-resolved. 
Velocity -unresolved observations toward bright FIR contin¬ 
uum sources with foreground diffuse clouds thus result in ap¬ 
parently weaker [C n] luminosities. However, the most impor¬ 
tant effect revealed by the velocity-resolved [Cn] emission 
observations is the decrease of L[Cii]/Lfir by large factors 
as one moves from the extended cloud component to the FIR 
continuum peak. Just like we have demonstrated in Orion (the 
L[Cii]/L F ir and the f=[N(C + )/N u ]/lAxlO~ 4 maps shown in 
Figures [8jd) and(9le) are similar), we believe that the decrease 
of L[C ii]/L F ir is"due to the increase in the total column den¬ 
sity throughout the cloud relative to the (surface) [C n] emit¬ 
ting column. 

In addition, loca l low-metallicity galaxi es such as the 
nearby LMC (e. g., |Lebouteiller et al. 2012), but also early 
galaxies at z~6 (Capa k et al.||201 ffi T show enhanced [Cn] 
emission relative to their FIR emission. These values are in¬ 
terpreted as the consequence of higher penetration of FUV 
photons due to the lower dust content (lower extinction). This 
agrees with our view of galaxies having clouds of larger C + 
column densities relative to their total dust column. 


15 In the most general case, /[cii]//FiR-/[Cii](^k)/[ fB A (T d )(l - e T ^)dA], 
with and T d only coupled inside clouds at high densities n>10 6 cm -3 . For 
the lcr range of T d and Td,i 60 values toward OMC 1, 7 [Cii]//fir can be ap¬ 
proximated by 7 [cii]/7fir~C • 7[Cn]/[7?i6o(7d)(l - e -7 ^ 160 )], with the C con¬ 
stant varying a factor of -2 depending on the exact T d and Td,i60 values. 


Thus, if T d does not vary much, the FIR intensity raises with increasing total 
column density through the cloud, as shown in Figure|TT[b). 

16 At the average distance of -4 kpc of these massive star-forming regions, 
the ~50"x50" mapped areas are equivalent to -lxl pc 2 . The typical spatial 
resolution toward these more distant SFR regions with Herschel is ~0.2 pc. 
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Figure 16. ( a ) Simplified geometry (not to scale) of OMC 1 (red, molecular gas and dust), the FUV-illuminated cloud face (blue, dominating the [C n] emission), 
dense ionized gas in the Hu region surrounding the Trapezium (green), and the atomic Veil (shaded, producing foreground Hi absorption). The observed [Cn] 
emission, /[cnp reaches the observer (to the left) almost unattenuated. At a given position, the observed ~ 160/mi dust emission is 5i6o(Fd)(l - ^ _Td ’ 160 ), where 
t d,i60 is the -160//m dust opacity along the entire line of sight. ( b ) Mixed C + and dust grains along the line of sight. In this case /[cii]=fo,[Cii](l - ^ _Td ’ 160 )/Td,i60 ? 
where /o,[Cii] is the unattenuated [Cn] emission, (c) [Cn] sources embedded behind large columns of dust. In this case /[cii] =fo,[Cii] e _TcU60 . 


5.4.2. Dust Attenuation 


The Trapezium cluster and associated H n region are only 
half enveloped by the molecular cloud. The [Cn] emitting 
region is thus facing us and is mostly unattenuated by fore¬ 
ground dust. Other massive star-forming regions or galax¬ 
ies may be embedded in large amounts of dust (buried star- 
bursts in galactic nuclei in particular). In such cases, the 
[Cii] emission itself will be attenuated by e _TcU60 (a dust 
screen), where Td,i 60 now represents the dust opacity between 
the observer and the [Cii] emitting layers (Figure 16]:). For 
V(H 2 )>10 24 cm -2 (rd, 160^1), the expected L[Cii]/L F ir values 
can be very low, lower than 5xl0 -5 (black dotted curve in the 
inset of Figure [13]). The prototypical example in the Milky 
Way is Sgr B2, the most massive starburst in the Galactic Cen¬ 
ter. Sgr B2 also shows a pronounced decrease of L[C ii]/Lfir 
as the dust opacity increases from the extended cloud to 
the ( FIR) optically thick cores and embedded H u regions 
(e.g., [Goicoechea et al.|2004[|Etxaluze et al.|2013j ). This rea¬ 
soning agrees with recent observations of local LIRGs in 
which the L[C ii]/Lfir d eficit is restricted to thei r nuclei and 
not their extended disks ( [Diaz-Santos et al.|2014|) . 

An intermediate case may exist in which the C 4 gas and the 
FIR emitting dust grains are well mixed through the entire line 
of sight (Figure [T6b ). In this case, the intrinsic [C nJ_emission 
will be attenuatecTby (1 - e _T(U 60 )/Td,i 60 (e.g., Thronson et al. 
|1990| ). This model provides less extreme corrections at large 
extinctions (magenta dashed curve in Figure [15]). Given the 
surface origin of the [Cii] emission, this model is however 
less realistic for OMC 1, but it may apply to low metallicity 
galaxies (lower A v /Nh and dust-to-gas mass ratios) in which 
larger columns of C + relative to the total dust column exist due 
to the much higher penetration of the FUV radiation field. 

Protostellar sources in Orion BN/KL are embedded in 
large column densities of material and thus very low 
L[C ii]/L F ir ratios can be expected. Using the He rschel /PACS 
spectroscopic data (Goicoeche a et al. 2015| ) we derive 
L[C ii]/Lfir-2.4x10 -S and Td,i60-2.7 toward the hot core and 
IRc sources, and L[C ii]/Lfir-7.7x10 -5 and Td,i60-l-2 toward 
the adjacent H 2 Peak 1 outflow region. Both positions are 
shown in the inset panel of Figure [15] In the frame of the sim¬ 
ple geometrical models discussed above, the L[C ii]/Lfir ratio 
toward H 2 Peak 1 is still compatible with th z foreground sur¬ 
face [C n] emission model. The L[C ii]/Lfir ratio toward the 


hot core is even lower, but this is likely a consequence of the 
higher FIR luminosity due to internal dust heating from the 
protostellar sources. In any case, the embedded [C n] emis¬ 
sion model provides extinction corrections that are too large. 


5.4.3. Inhomogeneities and FUV radiation penetration 

The struct ure of molecular clouds is not perfectly homoge 
neous (e.g., |Stutzki & Guesten|[l990| [Falgarone et al.||l991 * 


and FUV radiation from massive stars can penetrate to deeper 
cloud depths than in homogeneous clouds. Our represen¬ 
tative (homogeneous) PDR models for /2 h=2x 10 5 ~ 4 cm -3 , 
and dust grain properties appropriate to Orion, predict 
A/XC + )^(1.2-1.6 )x 10 18 cm -2 . These C + column densities are 
a factor of ^2 larger than the predictions of PDR models us¬ 
ing standard ISM grain properties (Section | 3 . 2 | ). However, 
they are still lower, by a factor of - 3 , than the N( C + ) columns 
inferred toward many positions of OMC 1 (see Section [4]). 
These higher values suggest a degree of inhomogeneity at 
- 25 " scales (- 0.05 pc). As noted by Stacey et al7| jl 993 [ ), 

ss the pror ~~ r 


an 


inhomogeneous structure facilitates the propagation of FUV 
photons far from the Trapezium and contributes to the in¬ 
crease in the fraction of the volume that is effectively pro¬ 
ducing [C n] (more surfaces along the line of sight). 


6 . SUMMARY AND CONCLUSIONS 

We have presented the first ~7.5 , xll.5 / (~0.9pcx 1.4pc) 
velocity-resolved map of the [Cn]158/im line toward 
OMC 1. In combination with FIR photometry and H41<r and 
CO J- 2-1 line maps, we obtained the following results: 


1. The main contribution to the [Cii] luminosity (-85 %) 
is from the extended, FUV-illuminated face of the 
cloud (Go>500, tzh>5x10 3 cm -3 ) and from dense 
PDRs (Go>10 4 , cm -3 ) at the interface between 

OMC 1 and the H n region surrounding the Trapezium 
cluster. In addition, -15 % of the [Cii] emission arises 
from a different gas component without CO counter¬ 
part. Part of this emission can be associated with fila¬ 
mentary structures that photoevaporate from the cloud 
and with depressions of the visible-light emission previ¬ 
ously seen in the images of the Orion nebula (the Dark 
Bay, the Northern Dark Lane and other components in 
the atomic Hi gas Veil). An additional but minor con¬ 
tribution can be associated with ionized gas. 
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2. The highest C + column density peaks are found to¬ 
ward the Trapezium, the Orion Bar, and specific PDRs 
to the east of BN/KL. These bright [C n] emission re¬ 
gions (7p > 150 K) follow a spherical shell structure 
(in projection on the sky) surrounding the back edge of 
the blister Hu region. Several of these OMC 1/Hn in¬ 
terfaces have nearly edge-on orientations. The ex¬ 
citation temperatures and line opacities derived from 
the [Cn] and [ 13 Cn] lines are T QX = 250 - 300 K and 
T"[cii] - 1.5-3. Hence, the [Cn] line toward these 
peaks is slightly optically thick, and traces high exci¬ 
tation dense PDR gas (7T>300 K and tzh^IO 5 cm -3 ). 

3. L[Cii]/Lfir decreases from the extended cloud compo¬ 
nent (~10 -2 -10 -3 ) to the cloud center (~10 -3 -10 -4 ). 
The latter values are reminiscent of the “[Cn] line 
deficit” seen toward local ULIRGs. The L[Cii]/Lfir 
variations can be explained by the approximately face- 
on geometry of the molecular cloud relative to the ex¬ 
ternally illuminating stars. This creates a [C n] emission 
layer at the FUV-illuminated surface of the cloud. The 
FIR emission, however, approximately traces the bulk 
of material along each line of sight. Indeed, L[C ii]/Lfir 
is spatially correlated with the dust column density in a 
way that agrees with a simple face-on slab model. Try¬ 
ing to link the local and the extragalactic emission, we 
conclude that the [C n] emitting column relative to the 
total dust column density along each line of sight de¬ 
termines the observed L[C ii]/L F ir variations, with the 
lowest ratios observed toward the highest column den¬ 
sity peaks. This also implies that the specific grain 
properties and dust-to-gas mass ratio are important pa¬ 
rameters in determining the L[C ii]/Lfir ratio. 

4. We show that the T F iR/M Ga s ratio follows the hydro¬ 
gen ionizing photons traced by the H41 a line, reaching 
L F iR/M G as>80 L q M~ l toward the H n region around the 
Trapezium. These high L F iR/M Gas values are interesting 
proxies of the ionization parameter U, and are similar 
to the threshold where luminous galaxy mergers start 
to show “[Cn] emission deficits”. Indeed, the surface 
densities we infer toward OMC 1 (£ Gas -2000 M©pc -2 
and X S fr=(2.3-3.4)x10 -5 M 0 yr -1 pc -2 ) are compara¬ 
ble to those inferred in galaxies hosting vigorous star- 
formation. The interesting point of course is that these 
conditions, those of a young massive cluster still sur¬ 
rounded by its dense parental molecular cloud, can 
dominate the FIR emission at kpc galactic scales. 


The [Cn] 158 jam line is undoubtedly a bright and power¬ 
ful diagnostic of the FUV-illuminated ISM. In combination 
with other tracers ([Oi], CO, dust SEDs, PAHs and Hi), a 
very detailed picture of the dominant environment and pre¬ 
vailing physical conditions can be extracted. In addition, ob¬ 
servations of molecular ions (CO + , HOC + , CF + or CH + ) di¬ 
rectly related to C + , are becoming interesting tools to con- 


H 2 gas (e.g., Storzer et al. 1995 

Fuente et al. 2003; Neufeld 

et al. 2006; Guzman et al. 2012; 

vfagy et al. 2013). Even more 


low H 2 fraction (nearly atomic) gas where [C n] is also ex¬ 
pected ([Neufeld et al.|2010[ [van der Tak et al.|2013[ |Schilke 
|et al.|2014|). Most of them can be observed with ALMA. 
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Notes Added after Proofs : Three relevant papers were dis¬ 
covered after proof reading, below we include their most per¬ 
tinent information related to this work: 

• O’Dell & Harris (2010, AJ, 140, 985) presented a 
more detailed geometrical model of the Orion molec¬ 
ular cloud, Orion nebula and Veil (see their Figure 13). 

• Diaz-Santos et al. (2013, ApJ, 774, 68) determined the 
L[Cii]/L F ir ratio in a sample of 202 local LIRG and 
ULIRG galaxies. The L[Cii]/Lfir ratio was found to 
scale with the strength of the the 9.7 jam silicate absorp¬ 
tion feature, except for “optically thick” galaxies show¬ 
ing very large 9.7 /am strengths. Unlike the MIR dust 
emission/absorption, the FIR/submm emission is more 
sensitive to the total column density of dust in dense 
clouds. Thus, their observed trend is consistent with 
our detailed interpretation of the L[C ii]/Lfir variations. 

• Compared to local ULIRGs, Rigopoulou et al. (2014, 
ApJL, 781, LI 5) found enhanced L[Cii]/Lfir lumi¬ 
nosity ratios in a sample of intermediate redshift 
(0.2<z<0.8) ULIRGs observed with Herschel/ SPIRE. 
This result is more similar to observations of high red- 
shift galaxies with ALMA, and can also be understood 
as variations of the [C 11 ] emitting column relative to the 
total dust column density toward each galaxy. 
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APPENDIX 

APPENDIX: SED FITTING EXAMPLES 



Figure 17. Appendix: Example of SED fits to Herschel/ PACS-SPIRE and JCMT/ SCUBA2 photometric data in Jy/pixel units. All convolved at a uniform angular 
resolution of 25". Selected positions are indicated by their offsets with respect to the map center in arcsec. 








